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Preface 


The  AFCRL  program  of  radar  measurements  supporting  the  SAMS  rain 
erosion  project  at  Wallops  Island,  Virginia  during  the  1970-73  season  was  conducted 
as  a  joint  team  effort  by  various  members  of  the  Weather  Radar  and  Convective 
Cloud  Physics  Branches  of  the  Meteorology  Laboratory,  AFCRL.  Contractual 
support  was  provided  by  the  Applied  Physics  Laboratory  of  Johns  Hopkins  Univer¬ 
sity.  The  program  was  directed  by  Dr.  Robert  M.  Cunningham,  of  the  Convective 
Cloud  Physics  Branch,  and  supervised  by  Dr.  Kenneth  R.  Hardy,  of  the  Weather 
Radar  Branch. 

The  AFCRL  and  contract  contributors  to  the  radar  measurement  program  are 
identified  in  the  following  list,  and  their  efforts  toward  the  accomplishment  of  the 
SAMS  objectives  are  acknowledged.  In  addition,  thanks  are  extended  to 
Mr.  Alfred  A.  Spatola  and  Dr.  Kenneth  Hardy  of  AFCRL  and  Mr.  Hugh  Church 
of  the  Sandia  Laboratories,  Albuquerque,  New  Mexico  for  their  review  and 
criticism  of  the  manuscript.  Their  comments  and  suggestions  were  helpful  and 
appreciated. 
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A  Summary  of  the  Radar  Equations  and 
Measurement  Techniques  Used  in  the 
SAMS  Rain  Erosion  Program  at 
Wallops  Island,  Virginia 

AFCRL/SAMS  Report  No.  1 

I.  IN  THOM C.TION 

The  general  objective  of  the  SAMS  rain  erosion  program  at  Wallops  Island, 
Virginia  is  the  determination  of  the  erosion  effects  on  the  nose-cones  of  missiles 
caused  by  hydrometeors  located  along  the  trajectory  paths.  Missiles  are  launched 
into  storm  situations  that  move  across  the  Wallops  area  during  the  late  Fall, 
Winter  and  early  Spring.  Erosion  data  are  obtained  from  the  missiles  by  in-flight 
telemetry  and  by  nose -cone  recovery  and  examination,  following  the  flight. 
Comparison  radar,  aircraft  and  surface  data  are  acquired  concerning  the  type, 
size-distribution  and  liquid -water-content  of  the  cloud  and  precipitation  particles 
that  existed  along  the  trajectories,  or  in  their  vicinity.  The  erosion  and  latter 
data  are  then  correlated  in  attempts  to  explain  the  observed  erosion  effects  and 
to  design  theoretical  models  describing  the  physical  processes  involved. 

This  report,  which  is  the  first  of  a  series  describing  the  contributions  of  the 
Air  Force  Cambridge  Research  Laboratories  (AFCRL)  to  the  SAMS  program,  is 
concerned  with  the  quantitative  radar  measurements.  The  report  consists  of  two 
major  pa  *ts.  The  first  part  describes  the  radar  concepts,  definitions  and 
equations  ti.at  establish  the  background  of  the  measurement  and  analytical 


(Received  for  publication  28  January  1974) 
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procedures  involved.  The  subject  matter,  although  primarily  based  on  radar 
theory  and  literature  references,  is  specifically  directed  toward  the  SAMS  program, 
in  the  sense  that  the  radar  systems  and  parameters  discussed  are  those  im¬ 
mediately  pertinent  to  the  program. 

The  second  major  part  of  the  report,  commencing  with  Section  .'1,  describes 
the  particulars  of  the  Wallops  Island  radars,  the  calibration  methods,  the  missile- 
trajectory  measurements  and  the  analyses  that  were  employed  during  the  three 
seasons  of  SAMS  operation  since  1970. 

Specific  data  for  individual  missile  flights  are  not  supplied  in  the  present 
report.  The  single  exception  is  the  listing  of  the  radar  calibration  constants.  The 
particular  data  for  the  separate  flights  will  be  presented  in  later  reports  of  the 
series. 

Because  of  report  completion  deadlines,  the  subject  matter  herein  is  treated 
in  a  summary  fashion,  rather  than  in  comprehensive  detail.  Also,  regarding  the 
various  references  cited,  these  are  merely  particular  examples  of  the  work 
efforts  of  certain  authors  in  the  field.  Thus,  credit  reference  is  neglected  for 
numerous  other  authors  who  have  contributed  importantly  to  the  theoretical- 
experimental  background. 

Various  symbols  used  in  this  report  depart  from  the  standards  recommended 
by  the  AMS*  Committee  on  Radar  Meteorology.  Such  departures,  although 
minimized,  were  required  because  the  SAMS  program  encompasses  disciplines 
other  than  radar  meteorology.  Additional  parameters  will  be  defined  in  later 
reports  and  there  is  a  need  for  consistent,  non- redundant  symbology  throughout 
the  totality  of  the  descriptive  reporting. 

The  numerical  system  of  equation  reference  is  used  herein.  Thus,  the 
equations,  in  most  cases,  are  referenced  by  number  only. 

Forthcoming  reports  of  the  series,  as  presently  visualized,  will: 

(1)  Provide  background  information  and  explain  how  the  precipitation  liquid- 
water-content  and  other  paramsters  of  erosion  interest  were  determined  from  the 
radar  data. 

(2)  Describe  the  aircraft  and  surface  measurements  that  provide  knowledge  of 
the  types  and  quantities  of  the  cloud  and  precipitation  hydrometeors  that  existed  in 
the  Wallops  storm. 

(3)  Discuss  the  correlation  analyses,  involving  the  aircraft  and  radar  data, 
which  are  designed  to  improve  the  accuracy  of  the  empirical  equations  relating  the 
radar-reflectivity-factor  (Z)  and  the  precipitation  liquid -water-content  (M). 

(4)  Summarize  the  storm  situations  of  each  of  the  SAMS  seasons  of  operations 
and  present  the  data  results,  in  profile  and  tabular  form,  l~r  the  trajectory  paths 
of  the  eleven  SAMS  storm  missiles  launched  in  the  1970-73  season. 

American  Meteorological  Society. 
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These  reports  will  be  submitted  as  analytical  and  literary  progress  permits, 
not  necessarily  in  the  order  indicated  above. 

2.  Hi: VI Eft  OF  THE  PERTINENT  RADAR  CONCEPTS,  DEFINITIONS  AND  E OVATIONS 
RELATING  TO  THE  SAMS  PROGRAM 

Radar  equations  will  be  developed  in  this  section  which  pertain  to  single, 
point-source  targets  and  to  multiple,  hydrometeor,  scattering  targets.  Pulse 
integration  and  video  integration  will  be  described  and  the  radar  equations  will 
be  modified  for  application  to  radar  data  acquired  by  video  integration. 

2.1  Gain  and  lleamwidlh  of  a  Parabolic  Antenna 
2. 1. 1  ISOTROPIC  RADIATOR 

The  power  density  in  space  at  any  given  distance  from  a  source  of  microwave 
energy  that  radiates  isotropically,  is 


where  Pt  is  the  radiated,  or  transmitted,  power  at  the  source,  r  is  the  range 
from  the  place  of  power  measurement  to  the  source,  and  Sj  is  the  power  per  unit 
area  normal  to  any  radial  line  drawn  outward  from  the  source. 

2. 1.  2  RADIATOR  WITH  GAIN 

If  the  microwave  energy  emanates  from  a  transmitting  antenna,  or  array, 
which  is  shaped  or  arranged  such  that  the  power  density  in  the  transmission 
direction  exceeds  that  of  an  isotropic  radiator,  then  the  maximum  power  per  uni; 
area,  or  maximum  power  density,  is 


max 


PtGt 

4*r2 


(2) 


where  G.  is  the  antenna  gaLi  for  transmission  . 

1  1 
The  general  definition  of  G^,  for  e<ample  see  Kerr  ,  is  that  it  is  the  ratio 

of  the  power  radiated  by  an  isotropic  antenna  that  is  required  to  produce  the 

same  field  strength  at  the  same  range  as  that  produced  by  the  directional  antenna 

in  its  direction  of  maximum  transmission.  Thus, 


1.  Kerr,  D.  E.  (1051)  Propagation  of  Short  Radio  Waves,  McGraw-Hill, 
New  York. 
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where  is  the  average  field  strength  in  the  direction  of  maximum  transmission, 
and  dti  is  a  unit  solid  angle  viewed  from  the  transmitter.  If  the  direction  of 
maximum  transmission  is  defined  as  the  polar  axis  of  a  spherical  coordinate 
system  with  its  origin  at  the  antenna,  the  angles  O'  and  are,  respectively,  the 
zenith  angle  from  the  polar  axis  and  the  angle  from  this  axis  in  the  plane  per¬ 
pendicular  to  the  zenith  plane.  The  function,  f(o',d'),  is  the  " pattern  function" 
of  the  particular  antenna  and  is,  in  general,  a  complex  function  of  the  angles. 

The  term  |f(0',d')|  describes  the  range-normalized,  polar  diagram  of  the 
antenna. 

If  Sj  and  So,  in  Eq.  (3),  are  set  equal,  using  the  general  definition  of  gain 
stated  previously,  then,  3ince  the  numerator  term,  ,  is  equal  to  4 jr,  the 

equation  may  be  rewritten  as 

Gt  ~f  i  4*  1 2 -  •  (4> 

This  defines  the  transmission  gain  of  an  antenna. 

2. 1.3  DEFINITIONS  OF  BEAMWIDTH  AM)  GAIN  FOR  A  PARABOLIC 
ANTENNA 

If  the  transmitting  antenna  is  a  parabolic  dish,  as  is  the  case  with  the 
Wallopc  Island  radars  used  for  quantitative  storm  measurements  in  the  SAMS 
program,  the  beam  pattern  that  emanates  from  the  antenna  consists  of  a  very 
sharp,  well-defined  primary  lobe  of  radiation  plus  various  secondary  or  side 
lobes.  All  of  these  lobes  are  symmetrical  about  the  polar  axis. 

The  power  density  of  the  radiation  is  a  maximum  along  the  polar  axis  and 
quickly  drops  off  in  value  with  increasing  angular  distance  from  the  axis  (as 
viewed  from  the  antenna).  If  a  cross  section  is  drawn  through  the  primary  lobe, 
normal  to  the  polar  axis,  at  any  given  range  which  is  large  compared  to  the 
antenna  dimensions,  the  isolines  of  constant  pow’er  density  on  such  a  spherical 
section  are  approximate  circles  concentric  about  the  polar  axis.  The  particular 
circle  along  which  the  power  density  is  half  the  maximum  defines  the  "half 
power  circle"  of  the  radiation.  The  conic  surface  formed  by  the  integral  of 
all  such  circles  at  all  ranges  has  been  conventionally  employed  to  define  the 
so-called  "beam"  of  the  antenna.  This  cone,  or  beam,  has  a  total  angular  width 
0,  from  one  side  of  the  cone  to  the  other,  across  any  plane  containing  the  polar 
axis.  The  angle  0  is  referred  to  as  the  "beanwidth". 
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It  Is  presumed,  in  this  definition,  that  the  actual  lobe  pattern  of  radiation 
can  be  approximated  by  a  pattern  that  has  uniform  intensity  out  to  the  conic 
boundaries.  With  such  a  presumption,  it  may  be  demonstrated— for  example  see 

9 

Austin  and  Williams”  — that  the  transmission  gain  for  a  perfectly  constructed  and 
illuminated  parabolic  antenna  is  related  to  the  beam  width  as 

G.  =  -^-  .  (5) 

0 


I’robert-Jones^  has  discussed  the  several  erroneous  assumptions  that  led  to 
this  formerly -used  relationship.  He  has  demonstrated  that  the  theoretical 
relationship  between  the  gain  and  beamwidth  of  a  parabolic  antenna  is  more 
correctly  given  by 


G. 


2.-2 

*  Ka 


where 


K“  «(Jf2</5>  . 

a 


(6) 


(7) 


In  this  latter  equation,  <*  is  the  fraction  of  the  total  power  from  the  feed  inter¬ 
cepted  by  the  antenna  and  |i  is  a  factor  specifying  the  uniformity  of  the  antenna 

illumination.  For  a  perfect  antenna,  I’robert -Jones  states  that  a  (I  =  1.0, 

2  _ 

while  f  ((1)  1.04,  giving  K  1.02.  He  also  cites  other  theoretical  work  to 

3 

show  that  K  ,  in  the  ideal  case,  has  a  value  close  to  unity, 
a 

Of  course  actual  antennas,  which  arc  used  operationally,  depart  from  the 
ideal.  In  such  a  case,  the  parameter  Ka  may  be  regarded  as  a  "calibration 
constant"  that  specifies  the  actual  performance  relative  to  the  theoretical. 

Certain  waveguide  and  radom?  losses  of  the  transmission  system  may  also  be 
included  in  K  ,  depending  on  the  particular  calibration  methods  used. 

For  the  primary  X'adars  employed  for  the  SAMS  measurement  program  a* 
Wallops  Island  (reference  Section  3  for  specific  identification),  it  may  be  stated 
that  the  Ka  values  range  from  0.  88  to  0.  95. 


2.  Austin,  F.  M. ,  ami  W  illiams,  K.  L,  (1951)  Comparison  of  Radar  Signal 
Intensity  with  Precipitation  Hate.  Technical  Report  No.  14,  W  eather  Radar 
Research  Project,  M.  I.  T. ,  Cambridge,  Massachusetts. 

3.  Probert -Jones,  J.R.  (1902)  The  radar  equation  in  meteorology. 

Quart.  J,  Roy,  Meteorol,  Soc.  88:485-95. 
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2. 1.4  GAIN  IN  TERMS  OF  EFFECTIVE  ANTENNA  AREA 


An  alternate  expression  for  the  gain  of  a  parabolic  antenna,  which  is 

4 

important  to  the  subsequent  derivations,  is  (see  Silver  ) 

4  it  A 

Gr~~7~‘ 


(8) 


where  X  is  the  wavelength  of  the  radiation  and  Ag  is  the  "effective  area"  of  the 
antenna.  The  effective  area  is  related  to  the  actual  cross-sectional  area  of  the 
parabolic  dish  (across  its  face,  normal  to  the  principal  axis)  as 


A  =  K  A  . 
e  op’ 


(9) 


where  Ap  is  the  actual  area  and  Kq  is  a  constant  of  correspondence.  From 
Eqs.  ( 6),  (8)  and  (9),  it  follows  that 


K  = 


,2  ..2 

»rX  K 

a 


o  .  .  J  ’ 
4A  0 


(10) 


The  values  of  this  constant,  for  the  primary  radar;,  of  the  SAMS  program 
mentioned  earlier,  range  from  0. 44  to  0,  GC.  This  is  stated  without  proof,  but  is 
easily  demonstrated  from  the  radar  performance  values  listed  in  Tables  4,  5 
and  G, 


2.2  Single  Taqcrts  itnil  Continuous  Javo  Transmission 

2.2. 1  INCIDENT  AND  REFLECTED  POWER  FOR  SINGLE  ISOTROPIC 
AND  NON-ISOTROPIC  TARGETS  WITHIN  THE  BEAM 

If  a  single  target  of  cross-sectional  area  A^  is  present  within  the  antenna 
beam,  the  target  intercepts  an  incident  amount  of  microwave  power 


p.ciAt 

o 

4;r  r“ 


(11) 


4,  Silver,  S.  (1951)  Microwave  Antenna  Theory  and  Design,  McGraw- 
Hill,  New  York.  ' 
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If  the  target  has  a  reflection  coefficient p^,  and  if  it  reradiates  or 
"scatters"  the  incident  power  isotropically,  the  power  density  of  the  reradiation 
at  any  particular  range  r'  from  the  target,  is 


or,  from  Eq.  (11) 


P,OtPtA 


r  .2  2  ,2  * 

(4*)  r  r' 


If  the  target  does  not  scatter  isotropically,  it  is  conventional  to  define  a 
"backscatter  cross  section",  a,  to  replace  the  geometric  cross  section  p{At. 

The  values  of  a  will  depart  from  those  of  p^At  depending  on  such  things  as  the 
target  dimensions  relative  to  the  wavelength  of  the  radiation,  the  aspect  angle 
and  curvature  properties  of  the  target  in  the  beam,  and  the  wave  polarization. 
However,  it  is  always  possible  in  principle  to  determine  the  a  value  for  any 
prescribed  target. 

With  a  used  to  replace  p^Aj,  Eq.  (11)  becomes 
p  Ca 

S.  - - (14) 

4^pt 

and  Eq.  (13)  is  written  as 

p  tGta 

S  -  - 2  2  2-  <15) 

r  (47r)^r  r*^ 

2.2.2  POW  ER  RECEIVED  FROM  A  SINGLE  TARGET  UNDER  CONDITIONS 
OF  CONTINUOUS  WAVE  TRANSMISSION 

The  equation  development  is  best  illustrated  at  this  point  by  first  describing 
the  power  received  from  a  single  target  under  conditions  of  continuous  wave  (CW  ) 
transmission.  We  may  then  proceed  logically  to  the  more  complicated  situation 
of  pulsed  radar. 


*The  reflection  coefficient  for  metal  is  1.0;  for  other  dielectrics,  such  as 
water,  it  is  less  than  1. 0. 
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As  demonstrated  in  the  previous  section,  the  power  density  of  the  back- 
scattered  radiation  from  a  single  target  is  governed  by  Eq.  (13).  Then,  assume 
that  this  backscattered  radiation  is  received  by  a  second  parabolic  antenna  which 
is  located  adjacent  to  the  transmitter  and  is  pointing  in  the  same  direction  (toward 
the  target)  such  that  the  polar  axes  of  the  two  antennas  '  e  closely  parallel.  Also, 
assume  that  there  is  appropriate  shielding  between  the  two  such  that  none  of  the 
transmitted  radiation  directly  enters  the  receiver. 

If  the  effective  area  of  the  receiving  antenna  is  Ae’,  the  power  intercepted 
and  received  by  the  antenna  will  be 

P  =  S  A  1  ,  (16) 

r  re 


where  Sr  is  the  power  density  of  the  backscattered  signal  from  the  target  as 
specified  by  Eq.  (13).  With  Eq.  (15)  substituted  in  Eq.  (16), 


ptGt° 


, .  ,2  2  ,2 

(4tt)  r  r* 


-v 


(17) 


The  gain  of  the  receiving  antenna  Gr  is  related  to  the  effective  area  of  the 
antenna  Ae'  by  an  equation  analogous  to  Eq.  (8);  that  is. 


G. 


4jtA  ' 
_ £_ 


(18) 


This  equation,  when  solved  for  Ae'  and  substituted  in  Eq.  (17),  yields 


PtGtGrX‘o 

..  .3  2  ,2 
(4ir)  r  r' 


(19) 


which  specifies  the  power  received  from  the  target  under  CW  conditions  of 
transmission. 


2.3  Kudur  Equation  for  a  Single  Target 

Let  us  now  consider  the  situation  in  which  the  transmitting  and  receiving 
antennas  of  the  previous  example  are  one  and  the  same,  such  that 

G  G.  ,  (20) 

r  t 
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with 


r'  r  . 


Equation  (19)  then  becomes 


2 

PtG“  A  a 

7 7  TT 

(4ir)  r 


(21) 


(22) 


which  is  the  "radar  equation"  for  a  single  target  of  backscattering  cross- 
section  o. 

This  equation  for  single  targets  may  oe  employed  in  radar  calibrations  in¬ 
volving  the  tracking  of  a  target  of  known  dimensions  and  backscattering  proper¬ 
ties,  such  as  a  metal-coated  sphere  suspended  by  a  balloon  or  parachute.  Eater 
reference  will  be  made  to  this  equation  in  this  context. 

2.  I  Kudur  Pnmmi‘trr>  and  Equation-  fur  Multi |il Ibdrotm-ti-or  Targets 

2.4.  1  RADAR  Pl'LSE  LENGTH  AND  PULSE  VOLUME 

A  radar  systematically  transmits  pulses  of  brief  time  duration  at  a  particular 
pulse  repetition  frequency  (prf).  Part  of  the  power  in  these  pulses,  as  they  move 
out  from  the  radar,  is  intercepted  by  any  target  present  within  the  beam  and  is 
scattered  back  toward  the  radar,  where  it  is  detected  by  the  receiver  during  the 
time  intervals  between  pulses. 

The  radar  determines  the  distance,  or  range,  to  the  given  target  by 
measuring  the  time  interval  between  the  transmission  time  of  the  outgoing  pulses, 
to  the  reception  time  of  the  backscatter  return  pulses.  Thus,  the  range  to  the 
target  is 


where  At  is  the  time  interval  measured  by  the  radar  and  v  is  the  velocity  of 

10  -1  P 

microwave  propagation  (equal  to  2.0  X  10  cm  sec  ).  The  factor  2  in  the 
denominator  accounts  for  the  fact  that  the  radar  waves  travel  a  two-w-ay  path, 
to  the  target  and  return.  Hence,  the  range  is  half  the  total  path  length. 

If  the  time  duration  of  the  transmitted  radar  pulse  is  h1,  then  the  physical 
length  of  the  traveling  wavetrain  that  moves  outward  along  the  beam,  in  the 
r  direction,  is 


h  =  vph*  . 


(24) 
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Both  of  these  parameters,  h',  measured  in  time,  and  h,  measured  in  length, 
are  customarily  referred  to  as  ti  e  "pulse  length"  of  the  radar.  Sometimes  they 
are  also  referred  to  as  the  "pulse  width". 

The  volume  of  the  traveling  wavetrain  of  the  radar  pulse,  as  it  moves  out¬ 
ward  within  the  confines  of  the  conical  beam,  is 


\VT 


T<r  tan  o)  h 


(25) 


where  r  tan  Q  is  the  diameter  of  the  beam  across  any  major  chord,  normal  to  the 
polar  axis,  which  is  measured  at  the  mid-range  of  the  pulse.  If  0  is  a  small 
angle,  a  fraction  of  a  degree  to  several  degrees,  the  approximation 

0  r  tan  0  ,  ( 2f>) 


with  n  specified  in  radians,  is  valid  without  sensible  error,  such  that  Eq.  (25) 
may  be  simplified  to 


VWT 


7r(rs)2h 

~~i - 


(27) 


Without  detailed  explanation,  it  may  be  stated  that  the  return  pulses,  which 
are  scattered  back  to  the  radar  from  any  ensemble  of  targets,  emanate  from 
within  a  volume  section  of  the  beam  that  has  a  length  equal  to  half  the  pulse 
length.  This  occurs,  primarily,  because  the  transmitted  pulse,  and  the  definition 
of  pulse  length,  involves  a  train  of  waves  that  are  moving  in  a  single  direction, 
outward  from  the  transmitter.  The  backscatter  pulse  from  the  illuminated 
targets,  on  the  other  hand,  involves  the  two-way  propagation  of  radio  waves, 
from  transmitter  to  targets  and  return.  This  causes  the  time  duration,  or  length, 
of  the  return  pulse  to  be  half  the  transmitted  pulse.  (The  length  of  the  illuminated 
volume  containing  targets,  in  other  words,  is  half  the  pulse  length. ) 

Thus,  with  reference  to  Eq.  (27),  the  pulse  volume  of  a  radar  (having  a 
parabolic  antenna)  is  defined  to  be 

V  ,  I«T  ,  rtrrtfh  _  (28) 

2.4.2  ECHO  RETURN  FROM  AN  ENSEMBLE  OF  SCATTERERS  AND  THE 
CONCEPT  OF  VOLUME  REFLECTIVITY 

If  an  ensemble  of  scattering  targets  is  present  within  the  pui.se  volume,  each 
target  of  which  has  a  particular  radar  cross-section  a,  the  summation  of  tiic  ...inf 
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backscatter  signal  within  the  volume  is 


B  =  I'o.  .  (29) 

The  volume  reflectivity,  p,  is  the  summed  backscatter  cross  section  per  unit 
volume.  Hence, 

0  B/V  ,  (30) 


The  dimensions  of  q  are  L 

2.4.3  BASIC  RADAR  EQUATION  FOR  IA  DROMETEOR  PARTICLES 


By  replacing  the  o  parameter  of  Eq.  (22),  which  pertains  to  a  single  target, 
with  the  analogous  parameter  for  an  ensemble  of  hydrometeors,  the  basic 

radar  equation  describing  the  power  received  from  hydrometeors  may  be  obtained. 

Specifically,  if  Eq.  (31)  is  solved  for  Ea.  and  substituted  in  Eq.  (22)  to  replace  u, 

i  1 

2  2  2 

Pt  A  G.Vhn 

pr  =  * - TT~’  (32) 

r  512  jr  “r 


This  specifies  the  amount  of  power  that  will  be  received  from  hydrometeors 
within  the  pulse  volume  which  have  a  volume  reflectivity  p. 

It  should  be  noted  that  the  received  power  from  hydrometeors  decreases  as 
the  inverse  square  of  the  radar  range.  In  contrast,  see  Eq,  (22),  the  received 
power  from  a  single  target  decreases  as  the  inverse  fourth  power  of  the  range. 

2.5  Olkrr  Factors  Involved  in  the  Itudar  Fquuliun» 

There  are  five  other  factors  that  enter  the  radar  equations  in  their  fully 
developed  form.  The  first  four  are  factors  which,  depending  on  the  radar  wave 
length  and  conditions  of  operation  and  whether  the  targets  are  point-source 
targets  or  multiple  hydrometeor  targets,  act  to  decrease  the  received  power 
relative  to  the  amounts  specified  by  Eqs,  (22)  or  (32).  The  last  factor  is  a 
units  conversion  factor  that  is  required  because  of  the  non-consistent  set  of  para¬ 
meter  units  used  operationally. 
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Each  of  these  factors  will  be  discussed  in  turn  in  this  Section.  However,  the 
radar  equations  will  not  be  rewritten  at  each  step.  They  will  only  be  rewritten 
after  all  five  factors  have  been  described.  The  work  will  be  summarized  at  this 
point  and  the  radar  equations  will  be  presented  in  their  fully  developed  forms. 

2.5.1  I’ ROBERT -JONES  FACTOR 
3 

Probert -Jones  has  demonstrated  that  the  beamwidth  of  a  parabolic  (or  other) 
antenna,  for  the  radar  case  of  two-way  transmission  with  the  pulse  volume  filled 
with  hvdrometeors,  differs  importantly  from  the  beamwidth  as  defined  convention¬ 
ally  for  one-way  transmission.  His  work  was  based  on  a  detailed  consideration 
of  the  pattern  function  integrals  of  the  antenna,  see  Eqs,  (3)  and  (4),  and  of  the 
shape  of  the  power  density  profile  of  the  primary  lobe  of  radiation  that  existed 
within  the  conical  boundaries  of  the  antenna  beam  of  normal  definition. 

He  found,  in  essence,  that  the  "effective  solid  angle",  of  the  conic  sur¬ 
face  bounding  the  pulse  volume  — from  which  the  backscatter  contribution  to  the 
received  power  was  approximately  99.  7  percent  of  the  total  (the  remainder  being 
contributed  by  sidelobes)— was  related  to  the  solid  angle,  fi,  of  conventional 
beamw  idth  definition,  as 


n 


e 


.  n 

*  2  ln2  • 


(33) 


This  means,  from  the  definition  of  a  solid  angle,  that  the  "effective  beam- 
width",  0e,  for  two-way  radar  transmission  with  hydrometeor  targets,  is  related 
to  the  beamwidth,  o,  for  one  way  transmission,  as 

fle  SB02  ,  (34) 

where  Sg,  which  will  be  referred  to  as  the  Probert -Jones  factor,  or  "beam  shape 
factor",  is 

SB  TTHT  °-  721  *  (35) 


2.5.2  PARTIALLY -FILLED  PULSE  VOLUME 

If  the  pulse  volume  is  not  completely  filled  with  hydrometeors  and  if  there 
is  any  way  of  recognizing  such  situation  a-priori*,  it  is  possible,  in  theory,  to 


*  The  reference  is  to  situations  that  occur  at  the  side  and  top  boundaries  of 
storm  clouds.  The  radar  pulse  volume  at  the  boundary  range,  which  intercepts 
the  boundary,  contains  hydrometeors  in  one  part  of  the  volume  but  none,  only 
clear  air,  in  the  remaining  part. 
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introduce  a  correction  factor,  0,  into  the  numerator  of  Eq.  (32),  to  account  for 
the  fractional  reduction  of  received  power  that  would  occur. 

In  practice,  however— particularly  with  high  performance,  high  resolution 
radars  having  a  very  narrow  beam  and  small  pulse  volume— this  factor  is  seldom, 
if  ever,  considered,  and  is,  implicitly,  assumed  to  have  the  value  unity. 

The  assumption  that  il  =1.0  was  used  in  all  of  the  radar  work  and  analyses 
of  the  SAMS  program  at  Wallops  Island. 

2.  5.  3  ATTENUATION  FACTOR 

Another  factor  that  acts  to  reduce  the  received  power  in  precipitation  and 
cloud  situations  is  attenuation.  This  is  a  wavelength  dependent  factor  that  is 
normally  of  negligible  importance  at  the  longer  wavelengths,  longer  than  S  band 
(~  10  cm),  but  becomes  increasingly  important  as  the  wavelength  is  decreased, 
from  S  to  C  band  (~  5.  7  cm),  through  X  band  (~  3.  2  cm)  to  K  band  (~  0.  9  cm). 

Attenuation  occurs  because  the  microwaves,  in  traveling  to  any  given  range 
location  in  the  beam,  are  forced  to  traverse  a  two-way  path  through  intervening 
hydrometeors,  which  absorb  or  otherwise  (by  v,i  rious  scattering  mechanisms) 
cause  the  microwaves  to  lose  a  certain  fraction  of  their  power  per  unit  path 
length.  The  fraction  depends  on  the  type  of  hydrometeors,  on  their  precipitation 
rate  (or  liquid  water  content,  in  the  case  of  clouds)  and  on  the  radar  wavelength, 
as  previously  noted. 

The  attenuation  factor,  A,  which  is  introduced  into  the  numerator  of  Eq.  (32) 
or  Eq.  (22)  whenever  attenuation  is  a  matter  of  concern,  is  specified  by 

n.  n9  n 

A  =  aj  1  x  a2  £  x  .  x  1  .  (30) 

where  a^,  .....  a^  are  the  fractional  or  decimal  attenuation  values,  per  unit 

path  length,  which  pertain  to  the  different  types  (and  quantities)  of  the  hydro¬ 
meteors  that  exist  along  the  radar  beam  at  ranges  smaller  than  the  particular 

"target  range".  The  n^,  n,, . n^  exponents  of  the  equation  are  the  number  of 

unit  path  lengths  contained  in  the  two-way  traverse  distances  through  each  of  the 
hydrometeor  regions. 

Attenuation  equations  and  values  are  commonly  expressed  in  decibel  form. 

In  such  a  case,  the  total  dB  attenuation  is 

10  log  A  =  Lj  x  10  log  aj  +  x  10  log  a,,  + . +  L.  x  log  a.  ,  (37) 


where  10  log  a^,  10  log  a,,, .  10  log  a.  are  the  "attenuation  rates"  for  the 

different  hydrometeors,  in  dB  per  unit  path  length,  and  1.^,  L^, . I,j  arc 

the  two-way  traverse  distances  through  the  specific  regions. 
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Particular  attenuation  equations  and  values  are  presented  in  Table  1,  for 
rain,  snow,  water  cloud  and  ice  cloud.  Four  radar  wavelengths  are  considered 
and  the  equations  and  values  are  listed  for  temperatures,  precipitation  rates  and 
liquid  water  contents  that  are  typical  of  the  ones  observed  in  the  winter  storm 
situations  at  Wallops  Island  during  the  1972  and  1973  seasons.  The  tabulations  for  the 
different  radar  wavelengths  are  included  because  of  the  possibility  that  C,  X  or  K 
band  radars  might  be  employed  in  the  SAMS  program  in  the  future. 

The  typical,  maximum  values  of  attenuation  that  would  be  anticipated  in  the 
"average"  winter  storm  .at  Wallops  Island  are  summarized  in  Table  2.  Typical 
thickness  values  are  listed  for  the  different  hydronv;teor  regions  and  for  the  storm 
in  total.  Also  listed  are  the  two-way  traverse  distances  through  the  regions  for 
an  assumed  radar  elevation  angle  of  20°.  This  angle  was  chosen  because  the 
slant  ranges  through  the  regions  and  storm  at  such  angle  are  about  the  maximum 
to  be  expected  under  the  SAMS  conditions  of  operation,  't  he  attenuation  values  for 
these  maximum  traverse  paths  are  given,  in  dll,  for  the  same  radar  wavelengths, 
precipitation  rates  and  liquid  water  contents  as  specified  in  Table  1. 

The  table  shows  that  attenuation  at  S  band,  even  for  the  maximum  values  listed, 
is  negligible  and  can  be  ignored.  Attenuation  at  C  band  is  indicated  to  be  generally 
negligible,  except  under  conditions  of  long  traverse  distance  (to  target)  or  of 
atypically  severe  storms.  Attenuation  at  X  band  is  a  factor  of  importance,  which 
should  be  considered  in  the  radar  equations  and  not  neglected.  Attenuation  at 
K  band  is  very  severe,  so  severe,  in  fact,  as  to  suggest  that  quantitative  measure¬ 
ments  in  this  wavelength  region  would  be  extremely  difficult. 

2,3.4  S’!  'T2M  l.OSS  FACTOR 

There  are  various  losses  in  any  radar  system  that  will  cause  the  power 
received  from  single  targets  or  hydrometeors  to  be  less  than  predicted  from  radar 
theory.  Also,  there  are  various  uncertainties  concerning  the  true  values  of  the 
various  parameters  that  enter  the  radar  equations.  For  example,  it  is  difficult, 
if  not  impossible,  to  state  that  the  transmitted  power,  or  the  pulse  length,  or  even 
the  wavelength,  the  gain  or  the  beamwidth,  had  definite,  specific  values  on  any 
given  day  or  period  of  radar  operation. 

However,  if  the  radar  is  calibrated  as  a  system  entity  — as  is  t he  case  with 
the  Wallops  Island  radars  used  for  quantitative  measurements,  as  will  be  described 
in  a  following  section— it  is  |X)sslble  to  collect  all  of  the  system  'osscs  and  un¬ 
certainties  into  a  particular  "loss  factor",  I.,  This  factor— if  used  in  the  numerator 
of  the  radar  Eqs.  (22)  or  (32)  in  coniunction  with  the  particular  P^,  h,  >,  and  0 
values,  which  are  the  theoretical  values  for  the  radar  system  or  arc  any  refer¬ 
ence  values  which  reflect  the  nominal,  usual  radar  performance  — provides  the 
means  of  specifying  the  actual  system  performance  relative  to  the  theoretical  or 
reference  standard. 
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Table  1.  Attenuation  Due  to  Precipitation  or  Cloud,  for  a  One-way  Path.  Extracted  from  Table  13.  5  of  the  Handbook 
of  Geophysics  (1960).  The  decibel  values  of  attenuation,  10  log  a,  are  given  lor  the  temperature,  precipitation  rate 
and  liquid  water  content  values  that  are  typical  of  the  ones  observed  in  the  Wallops  Island  storms  of  the  1972  and 
1973  seasons.  The  corresponding  decimal  values  of  "a",  are  also  listed.  P  is  precipitation  rate,  in  mm  hr-1, 

M  is  liquid  water  content,  in  gm  m"3.  The  P  value  for  snow  is  in  mm  hr”l  of  melted  water 


The  nominal  performance  characteristics  of  the  Wallops  Island  radars  are 
listed  in  Table  3.  For  purposes  of  the  present  report,  these  listed  values  are 
defined  to  be  the  ones  that  specify  the  reference  performance  of  the  radars. 

Thus  the  loss  factor  L,  as  defined  herein,  is  the  system  calibration  factor 
(established  for  the  given  day  or  period)  which— if  employed  in  the  radar 
Eqs.  (22)  or  (32)  in  conjunction  with  the  reference  values  of  P  ,  h,  X,  and  o, 
stated  in  the  tables— will  permit  the  numerical  evaluation  of  the  constant  of 
relationship  between  the  volume  reflectivity  n  and  the  received  power  Pr  (for 
the  particular  day  or  period). 

These  loss  factors  will  be  discvssed  at  greater  length  in  Sections  3.3  and  3.4, 
which  are  concerned  with  radar  calibrations. 

2.  5.  5  UNITS  CONVERSION  FACTOR 

It  has  bee  i  conventional  in  the  past,  or  operationally  convenient,  to  employ  a 
variety  of  units  for  the  different  parameters  entering  the  radar  equation  which  do 
not  comprise  a  consistent  set.  For  example,  it  is  conventional  to  express  Pf  in 
watt 3  or  milliwatts,  X  in  meters  or  centimeters,  r  in  statute  or  nautical  miles, 
h  ir  microseconds  (whereas  *.  is  actually  a  length  unit  as  it  enters  the  radar 
equation),  and  n  in  m  1  or  cm  The  use  of  such  inconsistent  units  requires 
that  a  "units  conversion  factor",  U,  be  introduced  into  the  radar  equations  to 
establish  numerical  correspondence  with  a  common  system  of  units.  Additionally, 
it  should  be  noted  that  different  units  conversion  factors  will  be  required, 
depending  on  which  form  of  the  radar  equation  is  used;  that  is,  "^oint  target  form" 
versus  multiple-scattering,  "i.ydrometeor  form,"  and  on  the  specific  parameter- 
units  chosen. 

It  might  be  contended  that,  normally,  the  problem  of  units  conversion  should 
be  left  to  the  user  of  any  given  equation.  However,  in  the  case  of  the  radar 
equation— with  the  great  diversity  of  units  employed  operationally  and  in  view  of 
the  fact  that  improper  conversion  can  cause  large  <  rrors,  if  unrecognized— it 
was  deemed  essential  to  define  and  employ  this  factor  herein. 

2.6  Summon  and  Presentation  of  the  i'ulR-lleteloped  liquations 

To  summarize  the  five  factors  discussed  in  the  preceding  sections,  we  will 
rewrite  Eqs.  (22)  (for  a  point-source  target)  and  (32)  (for  hydrometeor  targets) 
with  all  pertinent  factors  included.  The  assumptions  made  for  the  SAMS  radar 
measurements  will  then  be  reviewed  and  the  equations  presented  in  their  fully- 
developed  forms. 
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In  the  case  of  a  point-source  target,  the  fully-developed  equation  becomes,  from 
Eq.  (22), 

PtG.2A2aALU 
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where  U  is  the  particular  units  conversion  factor  for  the  different  parameter  units 
employed.  The  factors  Sg  and  <1  do  not  enter  this  equation,  since  they  pertain  only 
to  hydrometeor  targets. 

For  hydrometeor  targets,  Eq.  (32)  rewritten  becomes 
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where  V  is  another  units  conversion  factor,  different  than  U  . 

1  3 

It  was  mentioned  in  previous  sections  that  Sg  -  and  (/,  1.0  for  all 

SAMS  work,  and  that  the  attenuation  factor  A  for  typical  winter  storms  at  Wallops 
Island,  Virginia,  was  generally  negligible  at  S  and  C  band  but  was  important  at 
X  and  K  band. 

W  ith  these  values  inserted  in  Eq,  (39),  the  radar  equation  for  hydrometeor 
targets  in  its  fully-developed  form,  is 

P.A2G.202hrjALU 

_ A _ L  /  A  f\\ 


1024  In  2»2r2 


where  A  is  either  included  or  excluded,  depending  on  the  particular  radar  em¬ 
ployed.  This  equation,  in  its  essential  respects,  was  originally  derived  by 
Marshall  et  al.  3 

2,7  l‘ul*p  InU'itmtion  and  \  idea  Inle^rutiun 


2.  7. 1  PULSE  INTEGRATOR 

In  the  cloud  and  precipitation  situations  of  the  atmosphere,  the  radar  volume 
reflectivity  rj  and  the  associated  received  power  Pr  vary  widely  in  value,  in  an 
Incoherent  manner,  both  temporally,  from  pulse  to  pulse  of  the  radar  at  any  given 
range,  and  spatially,  between  range-adjacent  pulse  volumes  during  any  given  pulse. 
Pulse-lo-pulse  and  range  variations  of  5  to  10  dB  in  Pr  are  common  and  usual  with 
hydrometeor  targets. 

5.  Marshall,  J.S. ,  I^ngille,  R.C. ,  and  Palmer,  W.  M.K.  (1947)  Measure¬ 
ment  of  rainfall  by  radar,  J.  Meteorol.  4(No.  6):  186-192. 


These  incoherent  signal  variations  result  from  the  constant  "position 
reshuffling"  of  the  hydrometeor  particles  within  the  radar  beam.  This  reshuffling 
occurs  because  of  the  differential  fall  velocities  of  the  particles  of  different  size, 
also  because  of  atmospheric  turbulence.  It  causes  the  summed  backscatter  cross 
sections  to  vary  from  pulse  to  pulse  and,  additionally,  there  are  reinforcements 
■*  cancellations  of  the  wavetrains  of  microwave  radiation  emanating  from  the 
pulse  volumes,  which  also  contribute  to  the  incoherent  fluctuations. 

Atlas6  has  shown  that  the  time  required  for  the  scatterers  to  redistribute 
themselves  sufficiently,  by  turbulent  and  differential-fall -velocity  motions,  to  give 
non -autocor related,  independent  data  is 

r  0. 00171  A  ,  (41) 

where  >  is  in  centimeters  and  r  is  in  seconds.  He  also  notes  that  the  equation 

values  are  in  good  accord  vlth  the  times  reported  by  other  investigators,  such  as 

7  8  9 

Bartnoff  et  al,  Fleisher,  and  Stone  and  Fleisher. 

An  instrument  component  of  a  radar  system,  called  a  pulse  integrator,  was 
developed  by  Williams10  to  provide  stable,  average  values  of  received  power, 

Pr,  from  the  pulse-variable,  instantaneous  values  of  P  .  The  integrator  provided 
an  electronic  average  for  the  pulses  received  from  a  range  interval  (range  gate) 
equal  to  one-half  the  pulse  length.  From  initial  data  obtained  with  the  instrument, 
Williams  and  also  Austin  and  Richardson,  11  found  that  integration  times  of  one-half 
second  or  longer  were  required  to  obtain  a  good  statistical  average.  This  cor¬ 
responded  to  an  average  of  about  250  separate  pulses  of  the  radar.  The  pulse 
length  of  their  radar  was  1  psec;  the  prf  was  500,  as  best  judged  from  their  paper. 


(5.  Atlas,  D.  (1964)  Advances  in  radar  meteorology,  from  Advances  in 
Geophysics,  H);3 17-478,  Academic  Press,  New'  York. 

7.  Bartnoff,  S. ,  Atlas,  D. ,  and  Paulsen,  W.  II.  (1952)  Experimental  Statistics 
in  Cloud  and  Rain  Echoes.  Proceedings  3rd  Weather  Radar  Conference,  McGill 
University,  15-17  Sept  1952,  G1  to  G7, 

8.  Fleisher,  A.  (1953.  1954)  Information  Contained  in  Weather  Noise.  M.  I.  T. 
Department  of  Meteorology,  Weather  Radar  Research  Technical  Report  No.  22  A.B. 

9.  Stone,  M.  I..,  and  Fleisher,  A.  (1956)  The  Measurement  of  Weather  Noise. 
M.I.  T.  Dept,  of  Meteorol,  Weather  Radar  Res,  Report  No.  26. 

10.  Williams,  E.  L.  (1949)  The  Pulse  Integrator.  Part  A:  Description  of 
the  Instrument  and  its  Circuitry.  Technical  Report  No.  8A.  Weather  Radar 
Research  Project,  M.  I.  T. ,  Cambridge,  Massachusetts. 

11.  Austin,  P.  M. ,  and  Richardson,  C.  (1952)  A  Method  of  Measuring  Rainfall 
Over  an  Area  by  Radar.  Proceedings  3rd  Weather  Radar  Conference,  McGill 
University,  1 5-1 7  Sept  1952,  D-13toD-20. 
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The  pulse  Integrator,  following  these  initial  experiments,  was  quickly  applied 
to  the  hydrological  problem  of  determining  areal  rainfall  amounts  from  radar 
measurements,  Austin  and  Richardson'*  and  Farnsworth  and  Mueller***  described 
particular  versions  of  the  instrument  that  were  designed  for  this  purpose.  Con¬ 
siderable  development  work  along  these  lines  was  also  accomplished  during  the 
1950's  and  19fi0's,  which  led  to  further  refinements  of  the  Williams'  instrument. 

The  first  pulse  integrators  were  essentially  analogue  devices  which  provided 
a  temporal-spatial  average  of  the  square  of  the  voltage  amplitude  of  the  instan¬ 
taneous  received  signals.  The  signals  were  averaged  for  a  selectable  number 
of  individual  radar  pulses,  n^,  and  number  of  range  cell  elements,  n^.  Each  of 
the  latter  had  a  length,  in  the  r  direction,  equal  to  that  of  the  pulse  volume, 
which  is  h/2.  Thus  the  integrator,  in  essence,  provided  a  measure  of  the  average 
received  power 


n  n  1 


p  m 


(42) 


where  E  is  the  voltage  amplitude  of  the  received  signals,  I  is  the  circuit 
impedance  and  y  is  a  particular  constant  reflecting  the  differences  between  the 
electronic  average  and  a  true  average. 

Because  of  the  practical  difficulties  of  obtaining  Pr  values  for  a  large  number 

of  independent  samples  spanning  the  total  azimuth  angle,  range  and  altitude 

regions  of  observational  interest  to  radar  meteorologists,  and  also  because  of 

the  degradation  of  signal  resulting  from  integration  by  analog  methods,  Marshall 
13 

and  Hitchfield  suggested  that  Pr  could  be  determined  by  counting  the  fraction  of 
signals  exceeding  preset  threshold  levels.  This  count,  assuming  that  the 
probability  distribution  of  voltage  amplitude  is  known,  is  directly  related  to  the 
average.  Confidence  limits  can  also  be  specified  in  terms  of  the  ratio  of  the 
number  of  pulses  to  the  number  of  independent  samples,  which  tells  how  well  the 
electronic  average  corresponds  to  the  "true  average".  *  Marshall  and  Hitchfield 
presented  equations  and  graphs,  based  on  the  so-called  "weak  law'  of  large 


*  The  true  average  is  the  arithmetic  average  over  time  approaching  infinity. 

12,  Farnsworth,  G.W.,  and  Mueller,  K.A.  (1953)  Radar  Rainfall  Area 
Integrator.  Proceedings  Conference  on  Radio  Meteorology,  Art.  Dv-3, 

University  of  Texas,  Austin. 

13.  Marshall,  J.S.,  and  Hitchfield,  W.  (1953)  The  interpretation  of  the 

fluctuating  echo  for  randomlv  distributed  scatteres.  Part  I,  Can.  J.  Phys.  31: 
992-994.  -  ~ 
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numbers",  of  mathematical  statistics— for  example  see  Lindgren  and  McElrath^  — 
which  described  the  relationships  between  the  probability  distributions,  the  radar 
pulses  counted  and  the  independence  of  the  samples.  Atlas6  has  reviewed  the 
subject  and  extended  some  of  the  concepts. 

1  ri 

Various  types  of  pulse  integrators  exist,  for  example  see  Schaffner  ‘  and 
Groginsky,  that  incorporate  different  methods  of  averaging  and  data  presentation. 
These  include  (1)  analogue  techniques  of  averaging,  (2)  the  classification  of  signal 
counts  exceeding  specified  thresholds  of  signal  intensity  and  (3)  computer  methods 
of  averaging,  classification  and  processing.  In  general,  though,  the  measured 
average  received  power  is  related  to  the  true  average  received  power  as 


(43) 


where  P  is  the  true  average  and  y  is  an  "averaging  factor"  that  depends  on  the 
rt 

nature  of  the  electronic  averaging  performed  by  the  pulse  integrator,  when 
coupled  to  a  particular  radar  receiver  and  when  the  threshold  is  adjusted  in  a 

17 

prescribed  manner.  For  example,  Kodaira  has  shown  that  when  a  logarithmic 
receiver  is  used,  P  departs  probabilistically  from  the  true  average  dependent 
on  the  number  of  independent  data  samples  averaged  by  the  integrator.  When 
integration  is  performed  for  a  very  large  number  of  samples  approaching  infinity, 
~  0.  5(S2  Prj  (that  is,  y  0.  562,  or  -2.  5  dB)*.  For  single  independent  samples, 
on  the  other  hand,  50  percent  of  the  received  power  measurements  exceed  a 
threshold  level  that  lies  1.8  dll  below  the  true  average  such  that  Pr  0.  661  Pr 
(that  is,  y  :  0,601,  or  -1.8  dll).  This  illustrates  the  magnitude  of  the  y  values 
for  one  type  of  receiver  and  set  of  assumptions,**  Independent  samples  are 
discussed  at  greater  length  in  Appendix  A. 


Kodaira's  equations  were  actually  written  in  terms  of  mean  square  voltage 
amplitude.  However,  since  the  circuit  impedance,  sec  Kq.  (42),  is  a  constant  Tor 
any  given  circuit,  the  above  discussion  can  be  expressed  equally  well  in  terms  of 
received  power. 

**Referencc  Austin  and  Schaffner for  discussion  or  the  dependence  of  y  on  the 
spacing  of  the  input  quantization. 

14.  Lindgren,  B.W. ,  and  McElrath,  (J.W.  (1963)  Introduction  to  Probability 
and  Statistics.  The  Macmillan  Company,  New  York. 

15.  Schaffner,  M.  (1966)  A  Digital  Sweep  Integrator  for  Weather  Radar. 
Proceedings  12th  Conference  on  Radar  Meteorology,  17-20  October  1966,  Norman, 
Oklahoma. 

16.  Groginsky,  II.  L.  (1966)  Digital  Processing  of  the  Spectra  of  Pulse 
Doppler  Radar  Precipitation  Echoes.  12th  Conference  on  Radar  Meteorology, 

17-20  October  196G,  Norman,  Oklahoma. 

17.  Kodaira,  N.  (1960)  The  Characteristics  of  the  Averaged  Echo  Intensity 
Received  by  the  Logarithmic  I.  F.  Amplifier.  Proceedings  8th  Weather  Radar 
Conference,  Am.  Meteorol.  Soc. ,  Boston,  Massachusetts,  pp.  255-261. 

18.  Austin,  P.M.,  and  Schaffner,  M.R.  (1970)  Computations  and  Experi 
ments  Relevant  to  Digital  Processing  of  Weather  Radar  Echoes.  Proceedings 
14th  Weather  Radar  Conference,  17-20  November  1970,  Tucson,  Arizona. 
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When  a  pulse  integrator  is  used  as  part  of  a  radar  system,  the  radar  equation 
specifying  the  measured  average  power  received  from  hydrometeor  targets 
becomes  — from  Eqs.  (40)  and  (43),  assuming  that  the  details  of  integration  are  in¬ 
corporated  in  the  factor  •>  — 


yl^  2C’.“02hnAl.l 
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where  Tj  is  the  spatial-temporal  average  of  the  volume  reflectivity,  corresponding 
to  that  of  the  integrator.  More  conveniently,  Fq.  (44)  may  be  written  as 

1>  —4^  (45) 

r~ 

where 

i\A2a202hAyU 

c  1 -  (40) 
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is  a  constant  that  is  herein  defined  to  be  "the  radar  constant  of  system  reference 
performance".  It  is  seen  that  the  "averaging  factor",  discussed  above,  is 
included  in  this  constant. 

It  might  be  noted  that  the  radar  equation  for  a  single  target,  Kq.  (38),  remains 
unchanged,  in  the  case  of  pulse  integration.  The  single  target  gives  a  steady, 
coherent  echo  (relatively),  hencc"!^,  P  and  either  the  average  value  or  the 
instantaneous  can  be  employed  in  the  equation. 

2.7.2  AFC  IU.  VIDIX)  INTEGRATOR 

The  video  integrator  is  a  particular  modification  and  developmental  extension 

of  the  pulse  integrator.  This  instrument,  in  several  versions,  was  designed 

19 

by  the  Weather  Kadar  Branch  or  AIX'RI.,  for  example  see  Glover,  in  con¬ 
junction  with  the  Raytheon  C  orporation.  Two  different  models  of  the  instrument 
were  employed  in  the  SAMS  program  at  Wallops  Island,  Virginia  from  1970  to  the 
present. 

The  AFC'R I,  video  integrator,  in  essence,  provides  "  range  normalized"  values 
of  a  decibel  parameter. .  I,  which  is  related  to  the  integrated  volume  reflectivity 
as 


l  +  10  log  n  .  (47) 


19.  Glover,  K.  M.  (1972)  A  Precision  Digital  Video  Integrator.  Proceedings 
1 5th  Weather  Radar  Conference,  Am.  Metcorol.  Soc.,  Boston,  Massachusetts, 
p.  193-198. 
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wherej^  is  a  dB  constant  that  depends  on  the  calibration  constant  of  the  radar 
(discussed  in  Sections  3.2  and  3.3)  and  also  depends  on  the  particular  "sensitivity 
threshold"  that  is  specified  for  the  integrator  when  coupled  to  the  radar  receiver.  * 
If  the  radar  equation,  Eq.  (45),  is  solved  for  'n  and  substituted  in  Eq.  (47) 
and  written  in  decibel  form,  we  obtain  the  "operating  equation  of  the  video 
integrator", 


This  shows  that  the  integrator,  in  its  effective  operation  and  dependent  on  the 

particular  calibration  and  threshold  values  of  L  and  k,  determines  the  pulse 

integrated  values  of  P  ,  provides  electronic  multiplication  or  "range  normaliza- 
2  r  — 

tion"  by  r  ,  and  presents  the  ,  I  ,  values  based  on  the  relationship  of  Eq.  (48). 

This  is  an  oversimplification  of  the  actual  operational  details  of  the  instrument 
but  basically,  the  dB  values  of  the  integrated  signal  J_,  are  obtained  in  this  manner. 

The  AFCRL  video  integrator  also  has  the  design  capability  to  categorize  the 
storm-received  signals  in  classes  of  threshold  exceedence.  For  example, 
seven  selectable  classes  were  available  in  the  video  integrator  used  in  the 
1971-72  SAMS  season.  These  were  set  at  threshold  levels  some  3  to  7  dB  apart 
(usually,  although  1  dB  intervals  were  also  possible),  which  would  subdivide  the 
total  ^  span  of  the  storm  or  cloud  system  into  seven  different  categories  of 
signal  intensity.  Thus,  with  the  radar  in  operation  in  the  RHI  or  PPI  scanning 
modes,  the  signal  levels  within  the  storm  could  be  established  for  these  seven 
categories  for  all  ranges,  elevation  angles  and  azimuth  angles  of  operational 
interest. 

The  video  integrator  additionally  permits  the  signal  level  classes  to  be  dis¬ 
played  on  RHI  or  PPI  scopes,  in  different  "gray  shade  levels"  of  video  contrast 
(black,  dark  gray,  light  gray  and  white,  with  the  same  shading  sequence  being 
repeated  once  again  in  an  analytically  unambiguous  fashion).  Color  contrasts, 
with  television  presentation,  can  also  be  used  to  differentiate  the  various  classes 
of  signal  level.  The  Weather  Radar  Branch  of  AFCRL  is  currently  designing  and 
contracting  for  a  new  video  integrator  that  will  piovide  as  many  as  sixteen  classes 
of  color  contrast.  This  integrator  will  be  available  lor  the  1974  season  of  SAMS 
operations. 


♦It  will  be  conventional,  herein,  to  use  the  "bracket  underline"  symbol 
"  i_j"  to  specify  the  decibel  equivalent  of  any  given  parameter.  Thus,  with 
reference  to  Eq.  (47),  Jj  and  usj  arc,  respectively,  the  decibel  equivalents  of 
the  actual  signal  value  T  and  the  actual  calibration  constant  k. 
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To  summarize,  the  AFCRL  video  integrator  provides  recorded  and  real  time 
information  about  the  internal,  volume -reflectivity  structure  of  storms  and  cloud 
systems.  The  received  signals  are  averaged  for  a  selectable  number  of  radar 
pulses,  n  and  for  an  also  selectable  number  of  range  cell  elements,  n^.  The 

averages  are  determined  and  categorized  by  signal  level  using  the  basic, 

13 

threshold -exceedence  techniques  of  Marshall  and  Hitchfield.  The  real  time  data 
are  graphically  portrayed  on  RHI  or  PPI  scopes,  in  contoured  form  with  video 
shading.  The  recorded  data  are  processed  by  computer  and  used  in  the  quanti¬ 
tative  analyses  of  the  SAMS  program.  The  nature  of  the  display,  recording  and 
processing  during  the  different  seasons  of  operation  is  discussed  in  Section  4. 

It  is  pertinent  to  note  that  for  purposes  of  the  subsequent  equation  deriva¬ 
tions  and  manipulations  herein,  considerable  simplification  is  achieved  by  writing 
the  equations  of  the  video  integrator,  Eqs.  (47)  and  (48),  in  their  basic,  non-dB 
form.  Doing  this,  Eq.  (47)  becomes 

r  =  k“.  (49) 

and  Eq.  (48)  simplifies  to 


where  I  and  k  are  the  actual  values  of  the  integrated  signal  and  calibration 
constant,  rather  than  the  dB  equivalents. 

2.  7.  3  THRESHOLD  OR  SCALE  SETTING  OF  THE  VIDEO  INTEGRATOR 

When  the  AFCRL  or  other  video  integrator  is  used  as  an  instrument  compo¬ 
nent  of  a  radar  system  (coupled  to  the  output  of  the  radar  receiver),  it  is  necessary 
to  establish  or  set  the  integrator  during  the  process  of  coupling  such  that  a  parti¬ 
cular,  constant  relationship, 


exists  among  the  integration  signal  I  ,  the  received  power  Pr  and  the  square  of 
the  radar  range.  The  constant  c  will  be  referred  to  as  the  "coupling  constant". 

In  the  SAMS  program  ut  WuUops  Island,  the  coupling  constant  was  established 
by  determining  the  particular  integration  signal  I*  that  corresponded  to  the 
minimum  received  power  Prmin>  an^  was  "threshold  detectable"  at  a  specified 
"normalization  range",  rn<  Under  these  coupling  conditions,  Eq.  (51)  was  written 
as 
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c 


■= - ~~T 


min 


(52) 


and,  since  I  *  P,.  .  and  r  were  all  known,  the  value  of  c  was  computable. 

*  ’  1  nun  n  *  r  • 

For  example,  with  the  FPS-18  radar  at  Wallops  Island  in  the  1971-72  SAMS 

_  _q 

season,  the  threshold  received  power  was  set  at  -90  dBm  (Prnij„  '  10  m\V) 
at  r  -  1.0  NM.  The  corresponding  value  of  I*  was  "set  to  be"  495.45  (or 
Jj-  -  2fi.  95  dll).  Under  these  conditions,  c  4.95  X  10**  m\V  'em  *(NM) 

The  discussion  of  the  scale-set  conditions  for  the  video  integrator  will  not  be 
carried  past  this  point  of  illustration,  at  the  moment.  The  specific  conditions 
for  the  different  Wallops  radars  and  SAMS  seasons  will  be  described  later  in 
Section  3. 

2.7.4  RESOLUTION,  REPRESENTATIVENESS  AND  UNCERTAINTY 
CONSIDERATIONS 

When  radar  data  are  acquired  by  use  of  a  pulse  or  video  integrator,  various 
factors  must  be  considered  to  permit  the  quantitative  assessment  of  (1)  the 
temporal  and  spatial  resolution  of  the  measurements,  (2)  the  independence  and 
representativeness  of  ttic  hydrometeor  samples  included  in  the  integration,  and 
(3)  the  variance  of  the  measurement  values  abort  the  probable  "true  lues". 

The  discussion  of  these  factors  is  quite  involved  and  hence  is  presented  in 
Appendix  A.  Comments  about  the  possible  uncertainties  of  the  SAMS  radar 
measurements  are  included  in  Appendix  B. 

Certain  of  the  principal  facts  and  findings  of  the  first  of  these  appendices 
are  summarized  below: 

(1)  Pulse  or  video  integration  with  hydrometeor  targets  is  necessary  in 
order  to  obtain  stable  values  of  P  or  I,  and  to  permit  signal  level  categorization 
for  purposes  of  real-time  video  display  of  storm  echoes  and  for  photographic 
and/or  magnetic  tape  recording. 

(2)  Stabilization  of  the  integration  signals  may  be  defined  as  the  minimiza¬ 
tion  of  the  variance  of  the  P  or  I  values,  which  is  achieved  by  averaging  the 
instantaneous  received  signals  over  a  number  of  radar  pulses  and  over  multiple 
range-cell  elements. 

(3)  Excellent  signal  stabilization  may  be  achieved  by  integration,  but  the 
consequence  is  that  the  temporal  and  spatial  resolution  of  the  integrated  data  arc 
substantially  degraded,  relative  to  the  resolution  values  without  integration, 

(4)  The  degradation  of  the  temporal  resolution  is  of  no  practical  importance 
to  the  SAMS  program,  since  the  integration  times  normally  employed  are 
relatively  short,  on  the  order  of  a  few  tenths  of  seconds  (constant  for  all  radar 
ranges),  which  would  seem  eminently  suitable  for  SAMS  purposes. 
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(5)  The  degradation  of  the  spatial  resolufion  is  a  matter  of  SAMS  concern, 
however,  for  two  reasons:  first,  that  the  integration  volume  (particularly  its 
range  extent  if  multiple  range-cells  are  integrated)  is  large  compared  to  the 
spatial  resolution  of  the  telemetry  data  from  the  missile,  the  erosion  data;  and 
second,  that  the  integration  volume  is  also  large  relative  to  the  sample  volumes 
of  the  aircraft  measurements  (which  are  correlated  with  the  radar  data  as  will 
be  described  in  a  subsequent  report). 

(G)  The  standard  deviation  of  the  radar  signals  received  and  averaged  by  a 
pulse  or  video  integrator  is  specified  approximately  by  L'q.  (A2G)  of  Appendix  A. 
The  equation  shows  that  for  any  given  radar  having  a  particular  wavelength  V, 
beamwidth  0  and  pulse  repetition  frequency  (prf),  the  standard  deviation  of  the 
P  or  I  values  measured  with  the  integrator  will  decrease  (a)  with  the  number  of 
pulses  integrated,  that  is,  with  the  integration  period,  (b)  with  the  number  of 
range  cells  included  in  the  averaging,  and  (c)  with  the  scan  rate  of  the  antenna 
and  wind  speed,  within  certain  limitations. 

(7)  For  a  radar  having  a  selectable  prf,  the  standard  deviation  of  the  pulse 
or  video  integrated  signals  will  be  minimized  by  the  selection  of  the  largest  prf 
combined  with  a  long  integration  period  and  will  be  maximized  by  the  reverse 
combination. 

(8)  The  average  level  of  received  power  or  1  ,  witli  pulse  or  video  integra¬ 
tion,  remains  the  same  within  the  variance  limits  prescribed  by  the  equation 
cited  above,  for  any  given  value  of  the  radar  pulse  length.  Changes  in  the  average 
power  level  can  only  be  effected  by  changing  the  pulse  length;  the  level  increases 
with  h,  both  for  P  and  I. 


2.8  The  R«,fl«*rli\tt>  Kticlor 

Up  until  now,  we  have  not  specified  the  exact  nature  of  the  hydrometeor 
particles  which  are  the  scattering  targets.  This  was  deliberate,  since  such  con¬ 
sideration  involves  factors  other  than  purely  radar  factors. 

As  can  be  seen  from  the  radar  equation,  Eq.  (44),  a  weather  radar  provides 
measurements  of  the  volume  reflectivity,  of  the  scatterers  contained  within  the 
radar  pulse  volume.  Thus,  a  calibrated  radar  will  provide  true  values  of  q  (or 
n,  in  the  case  of  pulse  or  video  integration),  irrespective  of  the  nature  of  the 
scatterers.  The  scatterers  could  be  water  droplets,  ice  particles,  or  even 
insects  or  birds.  The  radar  in  each  case  will  give  an  accurate  measure  of  the 
summed,  backscattering,  cross-section  per  unit  volume,  which  is  q.  It  is 
essential  that  this  fact  be  clearly  understood.  To  reiterate,  the  fundamental 
atmospheric  parameter  measured  by  a  weather  radar  is  q. 

If  we  have  independent  knowledge  of  the  nature  of  the  scatters  being  observed 
by  a  radar,  microwave  scattering  theory  can  then  be  used  to  establish  the 


relationships  between  the  quantity  and  size  distribution  of  the  scatterers  and  the 
radar  volume  reflectivity.  In  the  case  of  atmospheric  clouds  and  storms,  the 
scattering  particles  will  consist  of  water  droplets,  for  temperatures  above 
freezing,  and  ice  particles,  for  temperatures  below  freezing.  There  is  the 
additional  case  of  supercooled  clouds,  which  can  exist  in  water  form  to  tempera¬ 
tures  well  below  freezing. 

For  Rayleigh  scattering  from  an  ensemble  of  water  droplets  (rain  or  water 

cloud),  assuming  no  coherent  scattering  and  ignoring  radar  pulse-shape  con- 
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siderations,  the  radar  volume  reflectivity  is  given  by  (see  Ryde  ) 


N(D.)D  j  , 


(53) 


where  N(Dj)  is  the  number  of  droplets  per  unit  volume  which  have  the  diameter 
Dj,  and  k  is  a  particle  shape  factor  that  considers  bulk  departures  of  the  water 
drops  from  spherical  shape.  The  factor  |  k|  ^  is  a  refractive  index  factor, 
specified  by 


(54) 


where  c  is  the  dielectric  constant  of  water  and 
m  =  n  -  jk'  , 

where  n  is  the  refractive  index  of  water,  j  -1  and  k'  is  the  absorption 
coefficient  for  water. 

A  radar  reflectivity  factor, 

Zw 

has  been  defined  by  Ryde  which  permits  Eq,  (53)  to  be  expressed  more  con¬ 
veniently  as 


n  r'^T"k  I Kl 2 Z\v  •  (57) 

20.  Ryde.  J.W.  (1946)  Echo  Intensities  and  Attenuation  Due  to  Clouds,  Rain, 
Hail.  Sand  and  Dust  Storms  at  Centimetre  Wavelengths,  GEC  Report  No.  7831, 
October  1941,  also  GEC  Report  No,  8516,  Aug.  3,  1944,  by  J.W.  Ryde  and 
D.  Ryde,  corrections. 
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The  shape  factor  k  in  this  equation  is  normally  assumed  to  be  unity,  that  is, 
the  water  drops  are  assumed  to  have  spherical  shape. 

The  refractive  index  factor,  |  k|  2,  for  water  (see  Gunn  and  East2^)  varies 
by  about  3  percent  with  the  radar  wavelength  over  the  X  range  from  X  band 
(~3.2  cm)  to  S  band  (~10  cm).  It  also  varies  by  about  0.5  percent  with  tempera¬ 
ture,  over  the  temperature  range  from  0°C  to  20°C. 

If  we  assume  an  average  temperature  of  10°C,  the  Gunn-East  data  give 
values  of  |  k|  2  of  0.  931  at  S  band,  0.  928  at  C  band,  0.  915  at  X  band,  and 
0. 873  at  K  band. 

Since  most  of  the  SAMS  radar  measurements  at  Wallops  Island  were  made 
at  S  band,  it  was  assumed  herein  that 

|  k|  2  =  0.  93,  (58) 


With  this  assumption,  and  with  the  assumption  about  the  shape  factor, 
Eq.  (57)  converts  to 


0.  93  t 


"W 


(59) 


where  the  subscript  W  on  Z  denotes  that  this  reflectivity  factor  pertains 

specifically  to  water  droplets,  in  the  form  of  rain  or  water  cloud. 

For  Rayleigh  scattering  from  ice  particles  (snow,  graupel,  ice  crystals, 
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etc.),  the  radar  volume  reflectivity  (see  Ryde  )  is  given  by 


n 


k|K|2H.N(Dr  )D.6 

1  M  M 


(GO) 


where  Dj  is  some  physical  measurement,  imprecisely  defined,  of  the  actual  size 

of  the  ice  particles,  and  N(Dj)  is  the  number  of  particles  of  this  size  per  unit 

volume.  The  other  parameters  are  as  previously  defined. 

Since  the  refractive  index  factor  divided  by  particle  density,  that  is,  K/p„ 

22  23  1 

is  a  constant,  from  the  theory  of  Debye,  Mason  has  defined  an  effective 
density -normalized -refractive -index -factor 

21.  Gunn,  K.  L.S,,  and  East,  T.D.R,  (1954)  The  microwave  properties  of 
precipitation  particles,  Quart.  J.  Roy.  Met.  Soc.  80:522.  (427-8). 

22.  Debye,  P.  (1929)  Polar  Molecules.  Chemical  Catalog  Co. ,  New  York, 
Section  8.  1  and  Chapter  V. 

23.  Mason,  B.J.  (1971)  The  Physics  of  Clouds.  Second  Edition,  Clarendon 
Press,  Oxford,  England. 
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I  k|pw\ 


where  is  the  density  of  water  (equal  to  unity)  and  pj  is  the  density  of  solid  ice. 
By  use  of  this  factor.  Mason  writes  Eq.  (60)  as 

5 

n  *  JrklK'PE .  N(D  )D6  ,  (62) 

X  1  ei  ei 


where  Dg  is  the  melted,  or  water  equivalent,  diameter  of  the  ice  particles.  He 
then  defines  a  reflectivity  factor,  analogous  to  that  for  water. 


v^i^vv 


such  that  Eq.  (62)  becomes 


n  =-<-k|K'|2z  . 

A  1 


Since  |  k|  2  =  0. 176  for  solid  ice,  and  since  pf  =  0.  917  gm  cm" 
adjusted  refractive-index-factor  of  Mason  is,  from  Eq.  (61) 


lensity- 


0, 176 
0.9172 


=  0.209  . 


This  value,  when  inserted  in  Eq.  (64),  yields 


n  --Qi-2p9£lkz 

n  a  K  • 

x 


This  gives  n  for  a  distribution  of  ice  particles  in  which  the  particle  sizes  are 
measured  in  terms  of  their  melted-equivalent  diameters.  The  particle  shape 
factor  k,  in  lieu  of  other  knowledge— and  since  it  is  probabilistically  likely  that 
irregularly -shaped  particles  will  be  randomly  oriented  within  the  radar  pulse 
volume— is  customarily  assumed  to  be  unity. 

The  q  equations  for  water  and  ice,  Eqs.  (59)  and  (66),  when  written  in  terms 
of  the  measurement  units  conventionally  employed,  become  for  water 


2l  ^  J  y  I  Q  —  1  2  y 


39 


with  Z^,  in  mm6m-3,  X  in  cm  and  n  in  cm”1,  and,  for  ice,  with  k  1.0, 

0  -a^xio-12Zl, 


<G8) 


where  the  units  are  the  same. 

The  review  of  the  radar  concepts,  definitions  and  equations  pertaining  to  the 
SAMS  measurement  program  at  Wallops  Island  is  terminated  at  this  point.  The 
applications  of  the  equations  to  the  program  objectives  will  be  described  in  the 
following  sections. 


:l.  Till-:  K\I)\H  \IK\>I  KKMKM  I'HtM.KXM  VI  »  M.MII'S  ISI.WII.  MIM.IM  \ 

The  general  nature  of  radar  measurements  and  end-product  data  obtained 
at  V  allops  Island,  Virginia  will  be  discussed  first  in  this  section.  The  different 
radars  and  their  performance  characteristics  will  then  be  described  and  the  two 
methods  of  radar  calibration  will  be  explained. 

.'I.  I  The  Ix-ni-ral  Nature  of  the  Hudur  Measurement-,  and  the  Knd-Pmduel  Data 

During  each  of  the  first  three  SAMS  seasons,  the  general  procedure  has  been 
to  launch  the  test  missiles  from  the  coastline,  from  "launch  pad  zero"  (see  map 
Figure  1)  in  the  14G°  azimuth  direction.  Thus,  the  trajectory  paths  of  the  missiles 
were  contained  in  the  vertical  plane  of  the  14G°  azimuth  angle  from  pad  zero. 

The  details  of  the  two-stage  rockets  and  payloads  have  been  described  by  Cole 
and  Robinett,  and  Cole,  Church,  and  Marshall. 

The  146°  launch  direction  was  chosen  because  the  primary  radars  on 
Wallops  Island,  used  for  the  SAMS  measurements,  are  located  at  the  Radar 
Atmospheric  Research  Facility  *  (RARF)  site,  which  lies  in  the  32(>°  azimuth 
direction  from  pad  zero  at  a  distance  of  1.  8  NM.  This  is  the  reciprocal  of  the 
launch  direction.  In  other  words,  radar  measurements  made  in  Rill  scan  at  an 
azimuth  angle  of  140°  were  coplanar  with  the  trajectory  paths  of  the  missiles. 

The  primary  radar  measurements  during  each  of  the  seasons  were  obtained 
from  an  S  band  radar  operating  in  the  Rill  mode.  The  end-product  data  desired 
from  the  measurements  were  the  values  of  the  radar  reflectivity  factor  along 
.he  course  of  the  missile  trajectories  at  the  time  of  the  missile  firing,  to  the 
best  time  correspondence  possible,  dependent  on  the  RIII  sweep  rate  and  timing 


*  This  facility  is  operated  by  the  Applied  Physic.,  Laboratory,  Johns 
Hopkins  University.  The  facility  is  also  called  the  JAFNA  (Joint  Air  Force 
NASA)  site. 
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of  the  radar  antennas.  (The  hydrometeor  parameters  were  determined  from 
these  reflectivity-factor  values,  as  will  be  described  in  the  second  report  of  the 
series.)  Other  radar  data  of  auxiliary  type  were  additionally  acquired,  as  will 
be  explained  later. 

The  values  of  the  radar  reflectivity  factor  were  established  in  the  following, 
step-by-step  fashion. 

The  values  of  the  video -integrated  signal,  I,  along  the  missile  trajectories 
were  first  determined  from  the  launch-time  Rill  photographs  (in  the  1970-71  and 
1971-72  seasons)  or  from  magnetic-tape  records  (in  the  1972-73  season).  Then, 
from  Eq.  (49)  inverted,  that  is, 

7T-— .  (69) 

K 

and  with  knowledge  of  the  radar  calibration  constant  *,  the  values  of  the  volume 
reflectivity  were  computed  for  each  of  the  trajectory  paths.  Finally,  with 
cognizance  of  the  height  of  the  melting  level  within  the  Wallops  storms 
determined  by  radar,  radiosonde  or  aircraft,  the  reflectivity  factor  Z^., 
for  the  "water  hydrometeors'  present  in  the  lower  portion  of  the  storms,  was 
computed  from  the  n  values  by  use  of  Eq.  (67)  inverted,  that  is. 


1.08 


Likewise,  the  reflectivity  .‘nctor  Zj  for  the  "ice  hydrometeors"  existent  in 
the  upper  portion  of  the  storms,  was  computed  from  Eq.  (68)  inverted,  or 


4.78  X* 


Across  the  melting  zone  itself,  the  values  of  the  reflectivity  factor  were  in¬ 
determinate.  This  zone  has  a  finite  altitude  extent,  because  the  snow  particles 
falling  from  aloft  into  warmer,  above-freezing  temperatures  melt  gradually  with 
fall  distance.  The  zone  is  composed  of  a  mixture  of  water-coated  snow’  particles 
and  liquid  droplets,  and  the  radar  backscattering  properties  of  the  mixture 
cannot  be  determined  without  gross,  questionable  assumptions. 


;;<Radar  information  about  the  melting  level  can  be  acquired  because  this 
level  commonly,  but  not  invariably,  shows  up  on  the  RHI  or  other  scope 
presentations  as  a  distinctive  "bright  band".  Reference  page  66  for  additional 
comments  about  the  melting  level  and  bright  band. 


Comparing  Eq.  (69)  with  Eqs.  (70)  and  (71),  it  is  readily  seen  that  the  end- 
product  values  of  Z^,  and  Zj  could  have  been  computed  directly  from  the  T  data. 
The  average  volume  reflectivity,  TJ,  is  merely  an  intermediate  parameter  in  the 
calculations  and  was  computed  only  because  of  its  fundamental  importance  in 
radar  meteorology.  In  other  words,  o  may  be  eliminated  between  Eqs.  (69)  and 
(70)  to  obtain 
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(72) 


and  between  Eqs.  (69)  and  (71)  to  obtain 

Zj  X  1012  .  (73) 

It  K 

These  equations  show  the  direct  relationship  between  the  radar  signal  values 
measured  by  the  video  integrator  and  the  desired  end-product  values  of  Z^.  and 
Zj.  The  only  parameter  of  these  equations  not  yet  described  is  the  radar  calibra¬ 
tion  constant  k.  We  will  explain  how  the  values  of  this  constant  w  ere  determined 
following  the  description  of  the  Wallops  Island  radars  and  their  particular  per¬ 
formance  characteristics. 

3.2  The  Radars  and  Their  Performance  Characteristics 

Although  casual  reference  has  been  made  to  the  several  radars  employed  in 
the  SAMS  program,  we  have  not  specifically  identified  the  radars  nor  stated  their 
locations  and  performance  characteristics.  W  ith  the  background  terminology  and 
equations  now  established,  we  may  proceed  to  these  descriptions. 

Five  radars  have  been  employed  in  the  SAMS  program  to  date,  for  various 
purposes.  *  Their  specific  locations  and  nominal  performance  characteristics  are 
given  in  Table  3.  Their  map  locations,  except  for  the  WSR-57  radar  at  l’atuxent 
River,  Maryland  are  shown  in  Figure  1. 

The  primary  radars  used  for  quantitative  measurements  w’ere  the  FPS-18, 

S  band  radar  (employed  in  all  seasons  to  date),  the  RARF  X  band  radar  (utilized 
in  the  1970-71  season)  and  the  Spandar,  S  band  radar  (used  in  the  1971-72  and 
1972-73  seasons).  **  All  of  these  radars  are  located  at  the  RARF,  JAFNA,  site 
on  W'allops  Island,  referenced  previously. 

*These  do  not  include  the  radars  used  for  missile-tracking  or  aircraft 
tracking. 

**The  FPS-18  radar  is  also  called  the  "monopulse  radar". 
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(le  Target  U,  MtOXtO-22  I.«»X10'22  I.4M  x  10*22  «.«0xl0'22  ».4»0xl0‘22 

Irometeor  Targets  U  8.736  XtO"7  8.738  X  10*7  0.736XIO*7  8.736X  10*7  8.736X  10"7 


Table  3  (Contd. )  Nominal  Performance  Characteristics  of  Reference  for  the  Radars  Employed  in  the  SAMS 
Program.  The  values  stated  are  the  values  normally  used  in  the  SAMS  program,  although  various  ranges  and 
choices  could  be  selected  concerning  certain  of  the  parameters 


Video  integrators  were  used  as  system  components  for  all  quantitative 
measurements  obtained  from  the  FPS-18  radar  and  also  from  the  Spandar  radar, 
in  the  1972-73  season.  Several  types  of  integrators  were  used  and  improved, 
season  by  season.  Differe.it  methods  of  data  recording  and  analysis  were  also 
used,  as  will  be  explained  in  Section  4. 

In  addition  to  these  primary,  data -gathering  radars,  there  were  two  other 
radars  that  were  utilized  in  a  qualitative  manner,  for  storm  monitoring  and 
launch-planning  purposes.  These  were  the  ASR-7  airport  surveillance  radar 
located  at  W  allops  Station  about  6  miles  NNE  of  the  RARF  site  (see  Figure  1), 
and  the  VSR -57  radar  located  at  Patuxent  River,  Maryland,  about  60  N'M  to  the 
MW  .  The  PPI  pictures  of  the  ASR-7  radar  were  transmitted  to  the  RARF  site 
and  presented  on  closed-circuit  television  for  real-time  viewing.  Facsimile 
reproductions  of  the  WSR-57  radar  PPI  photos  were  available  at  fhe  l'.  S.  National 
Weather  Service  facility  at  Wallops  Station  and  were  used  for  mission-planning 
purposes. 

The  information  for  the  ASR-7  and  VSR- 57  radars  is  presented  herein 
because  questions  about  the  relative  detection  capabilities  of  these  radars  and 
the  primary  radars  occasionally  arise  during  planning  operations.  Table  3 
provides  approximate  answers. 

The  reader  will  note  that  certain  of  the  parameter  values  in  Table  3  have 
been  estimated  or  left  unstated.  This  reflects  a  lack  of  calibration  information 
or  knowledge  about  the  particular  parameter. 

II.. 'I  Sphere  (.alihraliun  of  the  I'riimm  Kudar.x 

One  method  of  calibrating  the  primary  radars  was  the  so-called  "sphere 
calibration  method".  A  metal-coated  sphere,  or  krown  diameter,  radar  reflec¬ 
tivity  and  backscattering  cross  section  was  either  carried  aloft  by  a  balloon  or 
suspended  from  a  parachute  (dropped  from  a  balloon  or  aircraft)  and  was  tracked 
by  the  radar.  The  range  of  the  sphere  from  the  radar  was  known,  from  the  echo 
distance,  and  the  calibration  procedure  involved  the  monitoring  and  recording  of 
the  power  received  from  the  sphere  as  a  function  of  range. 

In  the  performance  of  such  calibration,  the  equation  for  a  point-source 
target,  Eq.  (38),  is  rewritten  as 


where  k^  is  the  "sphere  calibration  constant",  specified  by 

ks  kri.  ,  (75) 
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where  L  is  the  system  loss  factor,  as  defined  in  Section  2.  5.4,  and 


PtGt2*2AUs 

<4*)3 


<7«) 


is  the  constant  of  radar  reference  performance,  for  a  point-source  target. 

In  the  sphere  calibration,  the  values  of  kg  are  determined  directly  from  the 
measurements  of  received  power  versus  range.  Thus,  Eq.  (74)  is  written  as 


(77) 


and,  since  the  a  of  the  target -sphere  is  known,  and  since  paired  data  sets  of  the 
P  and  r  values  are  also  known  by  measurement,  the  particular  constant  value  of 
kg  can  be  determined.  It  might  be  mentioned  that  the  sphere  calibration  can  be 
performed  either  with  or  without  video  integration,  since  the  single,  spherical, 
target  scatters  coherently,  without  (appreciable)  pulse  to  pulse  variation  of 
received  power.  Hence,  Pr  *“  Pr,  and  either  can  be  employed  in  Eq.  (77). 

Once  kg  is  known,  it  is  a  simple  matter  to  ascertain  the  associated  values 
of  the  system  loss  factor  L.  The  constant  kr  is  first  evaluated,  from  Eq.  (76), 
for  the  reference  performance  values  listed  in  Table  3.  (In  fact,  kr  is  one  of  the 
listed  items. )  Then,  from  Eq.  (75),  written  as 

k 

L  =  t^-  ,  (78) 

r 

the  loss  factors  may  be  computed. 

Sphere  calibrations  were  performed  for  the  EPS- 18  radar  for  all  seasons, 
for  the  X  band  radar  in  the  1970-71  season  and  for  the  Spandar  radar  in  the  1971-72 
and  1972-73  seasons.  The  values  of  kg  and  L  for  these  radars  and  seasons  are 
shown  in  Table  4.  The  dates  listed  are  those  of  the  missile  launchings.  The  sphere 
calibrations  were  not  performed  on  these  exact  dates,  only  during  the  general 
periods.  However,  the  values  are  presumed  applicable  to  the  identified  days  of 
mlt  sile  firing. 

The  particular  calibration  constant  that  directly  enters  the  basic  computa¬ 
tional  equations  of  the  SAMS  radar  program,  see  Eqs.  (72)  and  (73),  is  the 
system  calibration  constant  k  .  The  values  of  this  constant  were  determined  from 
those  of  kg  as  follows. 
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From  Eq.  (50),  the  equation  of  the  video  integrator,  we  may  solve  for  k  to 
obtain 


C  LI 
_ r 

_  •)  ‘ 
P  r“ 
r 


(79) 


where  C  is  the  reference  constant  of  the  radar  equation  pertaining  to  hydrometeor 
targets,  as  given  by  Eq.  (40) 

Since  the  system  loss  factor  L  is  specified  by  Eq.  (78),  we  may  also  write 
Eq.  (79)  as 


K 


C  ck 
_E _ 3 


where 


(80) 


(81) 


is  the  "coupling  constant"  that  depends  on  the  particular  "scale  set  relationship" 
that  is  established  among  I,  Pr  and  r  when  the  video  integrator  is  coupled  to  the 
radar  receiver  and  used  to  obtain  measurements. 

The  AFC'RL  video  integrator,  as  discussed  in  Section  2.7.8,  was  coupled  to 
the  receiver  such  that  a  particular  integration  signal  7  was  received  in  associa¬ 
tion  with  a  specified  minimum  received  power  P r n , jj-, »  was  "threshold 

detectable"  at  a  normalization  range  rn<  1  NM).  Under  these  coupling  conditions 
the  constant  c,  of  Eq.  (81),  had  (was  "set"  to  have)  the  definite  value 


(82) 


Different  coupling  relations  were  used  for  the  different  radars  and  SAMS 
seasons  of  operations.  The  particular  Prmin>  rn«  T*  and  c  values  are  listed  in 
Table  4,  columns  6  through  9. 

The  desired  values  of  the  system  calibration  constant  k  were  computed  from 

the  k ,  values,  utilizing  Eq.  (80)  and  knowing  the  values  of  k  ,  C  and  c.  The 
s  x*  r* 

resultant  k  values  for  the  different  radars  and  seasons  are  presented  in  Table  4, 
column  11. 
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One  other  parameter  was  also  computed.  This  was  the  "calibrated  value"  of 
the  constant, 

C  =  CrL.  (83) 

of  the  basic  radar  equation  for  hydrometeor  targets,  Eq.  (45),  which  permitted  the 
specific  equations  for  the  individual  days  to  be  written  as 

P  (84) 

r 

The  C  values  are  shown  in  the  last  column  of  Table  4. 

3.  I  The  bi-ilromeler  Method  of  Radar  t'ali lira! ion 

The  second  method  of  radar  calibration,  which  will  be  called  the  "disdrometer 
method",  is  a  more  direct  and  straightforward  method,  equationwise,  than  that 
of  the  sphere  method,  but  it  is  also  considerably  more  difficult  to  accomplish, 
instrumentally  and  analytically. 

The  disdrometer  method  of  calibration  is  based  on  the  fact  that  the  system 
calibration  constant  k  can  be  computed  directly  from  Eq.  (72),  when  it  is  written 

as 


T 


where 

4  12 

1.08  \  X  10 

"  - 5 - 

it 


(85) 


(86) 


The  values  of  Z^.,  of  Eq.  (85),  were  determined  from  surface  measurements 
of  the  size  distribution  of  raindrops.  The  measurements  were  made  with 
"disdrometer  instruments",  as  will  be  described.  Simultaneous  radar  measure¬ 
ments  of  the  video-integrated  signal  T ,  in  the  numerator  of  Eq.  (85),  were  ob¬ 
tained  from  the  spatial  volume  of  the  storm  located  as  close  as  possible,  above 
the  site  location  of  the  disdrometers.  Independent  sets  of  Z^.  and  I  values  were 
acquired  for  a  period  of  dual  measurement  which  extended  an  hour  or  so  (usually) 
on  either  side  of  the  launch  time(s)  of  the  missiles.  These  sets  of  values  were 
then  employed,  in  conjunction  with  Eq.  (85),  in  comparison  and  correlation 
analyses  that  led  to  the  establishment  of  the  average  value  of  k  that  pertained  to 
the  calibration  periods. 
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Two  disdrometer  instruments  were  employed  for  radar  calibration.  (Other 
instruments  were  also  operated  but  for  other  purposes. )  These  two  disdrometers 
were  located  at  the  launch  site  of  the  missiles,  at  "pad  zero"  (see  the  map  of 
Figure  1).  The  second  instrument  provided  a  "data  consistency  check  and  back-up 
capability"  for  the  first.  These  disdrometers,  on  days  of  planned  missile  launching, 
were  operated  from  about  1  1/2  hours  prior  to  the  scheduled  launch  time  to  about 
112  hours  following  launch. 

Comparison  radar  data  were  obtained  throughout  the  operational  period  of  the 
disdrometers.  The  radar  being  calibrated  was  periodically  aimed  in  the  146° 
azimuth  direction  of  the  disdrometer  site  and  Rill,  A  scope,  or  other  data  were 
obtained.  From  analyses  of  these  data,  the  I  values  were  determined  which 
prevailed  about  1.3  km  offshore  of  the  disdrometer  site  at  altitudes  about  0.2  to 
0.  4  km  above  the  sea  surface.  Radar  data  could  not  1  ■  .  btained  from  directly  above 
the  disdrometer  site,  because  the  coastline,  the  launch  towers  and  the  other 
buildings  at  the  site  produced  "ground  clutter  return"  that  masked  the  desired 
precipitation  echoes.  Radar  data  had  to  be  obtained  from  beyond  the  shoreline  and 
the  1.3  km  distance -discrepancy  between  the  site  locations  of  the  disdrometer  and 
radar  measurements  had  to  be  tolerated  in  the  calibration  . 

With  these  introductory  comments,  we  may  now  describe  how  size-distribution 
and  Z„.  information  was  acquired  from  the  disdrometer  instruments. 

The  disdrometers  used  in  the  SAMS  program  were  designed  and  constructed  by 
24-26 

Joss  and  Waldvogel.  These  instruments,  in  essence,  measure  the  impact 

momentum  of  raindrops  of  the  various  sizes  that  fall  on  an  exposed  sensing  head 
2 

(of  50  cm  area).  *  The  instrument  is  primarily  sensitive  to  the  vertical  component 
of  momentum  resulting  from  the  gravitational  fall  or  the  drops,  and  is  relatively 
insensitive  to  the  horizontal  components  of  momentum  associated  with  wind  and 
turbulent  transport. 

-The  sensing  head  consists  basically  of  a  receiving  body  of  inverted,  conical 
shape  which  is  composed  of  styrofoam  and  capped  by  an  aluminum  disk.  This 
receiving  body  is  spring-suspended  within  a  containing  cylinder  such  that  the 
capping  disk  is  horizontally-exposed  to  the  precipitation.  Two  coils  are  attached 
to  the  lower  part  of  the  receiving  body,  their  axes  vertical.  These  coils,  on 
droplet  impact,  are  forced  downward  through  a  surrounding  electromagnetic  field. 
The  current  induced  in  the  lower  coil  is  amplified  and  used  to  provide  a  counter 
restoring-current  and  force  in  the  upper  coil.  The  amplitude  of  the  current  in  the 
upper  coil  is  a  measure  of  the  drop  size. 

24.  Joss,  J,,  and  W aldvogel,  A.  (1970a)  A  Disdrometer  for  Raindrops. 
Instruction  Manual,  Marc  Weibel  Dipl.  Ing. ,  ETH,  Kapellenstrasse  20, 

4000  Basel,  Switzerland. 

25.  Joss,  J,,  and  Waldvogel,  A.  (1970b)  Disdrometer  RD-69.  Instruction 
Manual.  Marc  Weibel  Dipl.  Ing.,  ETII,  Kapellenstrasse  20,  4000  Basel, 
Switzerland. 

26.  Joss,  J.,  and  Waldvogel,  A.  (1970c)  Analyses  AD-69  Instruction  Manual. 
Marc  Weibel  Dipl.  Ing,,  ETH,  Kapellenstrasse  20,  4000  Basel,  Switzerland. 
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The  electronic  measure  of  impact  momentum,  assuming  careful  calibration  of 
the  instrument  components  (sensing  head,  processor  and  recorder),  is  converted 
into  counts  of  the  number  of  drops  in  20  diameter  classes  (as  many  as  20,  depending 
on  the  size  range  of  the  rain)  that  have  fallen  onto  the  sensing  head  in  a  fixed  time 
interval.  To  be  more  specific,  the  instrument  provides  a  measure  of  the  number 
flux  of  the  drops,  Nj(Dj),  of  a  classified,  mid -diameter-size  D.,  that  have  fallen 
on  the  sensing  area  Ag  in  the  time  interval  Atg. 

The  rainfall  rate,  R,  is  directly  determinable  from  these  number  flux  measure¬ 
ments.  It  is  given  by 


R 


i  =  D 


— —  £ 
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i  =  D 


mm 
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(87) 


min 


where  D.  is  measured  in  mm,  D  .  and  D  are  the  minimum  and  maximum 
l  ’  mm  max  2 

diameters  observed,  and  Ag  is  specified  in  cm  with  Atg  in  seconds. 

The  number-flux  datr?  acquired  by  disdrometers  may  be  utilized  to  reconstruct, 

or  synthesize,  the  probable  concentration-and-size-distribution-properties  of  the 

raindrops  in  space,  as  they  existed  per  unit  volume  above  Hie  uisdrometer  sensing 

5 

head  just  prior  to  their  gravitational  fall  onto  the  head.  Marshall  et  al  have 
described  the  mathematics  of  this  reconstruction,  which  involves  an  assumption 

about  the  fall  velocity  of  the  raindrops  as  a  function  of  their  diameter. 

97  28  29 

The  experiments  of  Spilhaus,  Gunn  and  Kinzer,  I.iu  and  Orville  and 

others  have  shown  that  the  fall  velocity  of  raindrops,  within  the  restricted  size  range 

of  normal  drop  sizes,  varies  with  diameter  approximately  as 


,,  ,,b  cm 

VR  =aDi 


(88) 


for  D^  (classified  data)  specified  in  mm.  The  constant  a,  depending  on  investigator, 
varies  from  about  300  to  500  in  equation  consistent  units,  and  the  exponent  b  varies 
from  about  0.  5  to  0.  8, 


27.  Spilhaus,  A.F.  (1948)  Raindrop  size,  shape,  and  falling  speed, 

J.  Meteorol.  5:108-110, 

28.  Gunn,  R. ,  and  Kinzer,  G.  D.  (1949)  The  terminal  velocity  of  fall  for 
water  droplets  in  stagnant  air,  J.  Meteorol.  6:243,  (565,  594,  596-7). 

29.  I.iu,  J.Y.,  and  Orville,  If. D.  (1968)  Numerical  Modelling  of  Precipitation 
Effects  on  Cumulus  Cloud.  Report  68-9  Inst.  Atmos.  Sci. ,  South  Dakota  School  of 
Mines  and  Technology,  Rapid  City,  South  Dakota. 


From  this  equation  for  fall  velocity,  the  volumetric  concentration  of  the 
raindrops  in  space,  N(D.),  can  be  deduced  from  the  number  flux  data  provided  by 
the  disdrometer  instruments.  Thus,  for  any  given  diameter  class  D., 


N.(D.) 

N(D  )  =  — - — - —  -  -ff°» 

A  At  V_  T 
s  s  R  cm 


or,  from  Eq.  (88), 
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where  N(D.)  in  the  latter  equation  is  expressed  in  the  conventional  units  of  number 

3  1 

per  m  . 

With  knowledge  of  \(Dj)  thus  determined,  the  values  of  the  radar  reflectivity 
factor  for  water  drops  could  be  computed  by  writing  Eq.  (58)  as 


max  fi 
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and  employing  Eq.  (90)  to  obtain  the  computational  equation 
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The  subscript  "D"  on  Z^.  signifies  that  this  is  a  particular  reflectivity  factor, 
computed  from  disdrometer  measurements. 

The  Z\Vq  values,  in  practice,  were  determined  by  computer  from  magnetic  - 
tape  records  of  the  number  flux  data  of  the  disdrometers.  Numerous  other  cloud 
physics  and  radar  parameters,  such  as  the  precipitation  liquid  water  content 
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where  p^,  is  the  density  of  water  in  gm  cm  \  were  also  computed,  as  will  be 
mentioned  in  subsequent  reports. 

With  this  background  about  the  disdrometer  instruments  and  about  the  Z\vD 
computations,  we  may  now  return  to  the  subject  of  radar  calibration. 

It  should  first  be  noted  that  the  only  RARF  radar  that  was  calibrated  by  the 
disdrometer  method  was  the  FPS-18.  Calibrations  were  performed  for  the 
1970-71  and  1971-72  seasons  but  none  was  accomplished  for  the  1972-73  season. 
Initial  attempts  were  made  to  calibrate  the  FPS-18  radar  for  the  1972-73  season 
but  the  results,  at  least  for  the  day  of  27  February  1973,  revealed  an  apparent, 
substantial  discrepancy  between  the  disdrometer  and  sphere  method  results. 

This  discrepancy  has  not  yet  been  resolved.  It  is  suspected  that  the  disdrometer 
data  for  the  day  may  have  been  inaccurate,  due  to  lack  of  instrument  sensitivity  or 
possibly  because  of  improper  instrument  exposure,  (The  sites  and  exposure  were 
changed  between  the  1971-72  and  1972-73  seasons.) 

The  following  discussion  pertains  strictly  to  the  disdrometer  calibrations 
that  were  performed  for  the  FPS-18  radar  during  the  1971-72  season.  This  will 
serve  to  illustrate  the  general  method  of  calibration. 

The  radar  data  provided  for  calibration  purposes  in  the  1971-72  season  con¬ 
sisted  of  A  scope  values  of  the  video-integrated  signal,^  (in  dB)  that  were  ob¬ 
tained  from  the  "low  elevation  angle"  portion  of  RHI  sweeps  that  had  been  made  in 
the  146°  azimuth  direction  on  the  days  of  missile  firing.  The  data  were  read  from 
the  reccrded  A  scope  traces  at  radar  ranges  of  4.6  and  5.6  km  and  at  elevation 
angles  of  0.  5°  and  0.  75°.  The  four  values  obtained  were  averaged  and  this 
average  was  presumed  representative  of  the  radar  echo  conditions  in  the  rain  at  a 
distance  1.2  to  2.  3  km  offshore  of  the  disdrometer  site  at  altitudes  0.04  to  0.07  km 
above  the  sea  surface. 

These  average  Jj  values  were  obtained  for  most  Hill  scans  that  were  made  in 
the  146°  azimuth  direction.  The  values,  on  most  days,  were  acquired  for  the 
period  from  approximately  1  1/2  hours  preceding  launch  to  1  1/2  hours  after 
launch.  The  time  period  between  RHI  scans  varied  from  one  minute  (minimum)  to 
as  much  as  10  to  30  minutes.  This  depended  on  the  scanning  program  of  the 
particular  day.  * 

It  was  assumed  that  each  of  the  values  was  representative  of  a  one  minute 
time  period,  which  corresponded  to  the  minimum  spacing  period  between  the  radar 
data  points.  The  comparison  2\\  values  from  the  disdrometers  were  likewise 
computed  for  one  minute  periods  of  sampling. 


■'The  RHI  scans  with  the  FPS-18  radar  were  made  preferentially  in  the  146° 
azimuth  direction.  However,  during  the  aircraft  measurement  period  following 
the  missile  firing,  RHI  scans  were  also  made  in  the  azimuth  directions  of  the 
aircraft  locations,  which  sometimes  differed  from  146°. 
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Considerable  effort  was  expended  in  the  correlation  analyses  of  the  ^  and 
data  which  will  not  be  reported  in  detail  herein.  In  summary,  though,  there 
were  two  important  findings.  First,  it  was  discovered  that  the  correlations  and 
correlation  coefficients  of  the  data  sets  for  the  different  days,  were  appreciably 
improved  by  the  exclusion  of  the  radar  and  disdrometer  data  for  precipitation  rates 
smaller  than  about  0.  5  mm/hr.  The  comparative  data  for  rates  smaller  than  this 
evidenced  large  scatter  and  were  deemed  non-i  epresentative,  and  i  ere 

eliminated.  Second,  it  was  also  apparent,  because  of  the  distance  srj>ar.  :ion  of 
the  sites  of  the  disdrometer  and  radar  measurements  (1.3  km),  and  because  of 
the  occasional-to-frequent  presence  of  spatial  inhomogeneities  in  the  storm  rain¬ 
fall  at  this  distance  scale,  that  some  selection-rejection  criteria  were  required  to 
minimize  the  inhomogeneity  effects.  The  criteria  adopted  were  based  on  the 
assumption  that  temporal  uniformity  of  the  rainfall,  as  observed  at  the  disdrometer 
site,  was  likely  to  be  associated  with  spatial  uniformity  of  the  rainfall,  extending 
seaward  of  the  disdrometer  site  out  to  the  place  of  the  radar  measurements.  With 
this  assumption,  particular  data  points  of  J_,  and  Z\\D  were  eliminated  from  the 
data  sets  whenever  the  rainfall  rates  at  the  disdrometer  site  evidenced  large 
temporal  changes  that  departed  from  the  normal,  typical  variability  that  character¬ 
ized  the  particular  day. 

The  details  of  the  disdrometer  method  of  calibration  can  best  be  illustrated  by 
reference  to  a  particular  case  example.  The  day  selected  for  illustration  is  that 
of  17  February  1972.  The  dual  sampling  period  on  this  day  was  sufficiently  short 
to  allow  convenient  tabulation  of  the  data.  The  day  was  also  one,  incidentally, 
when  two  separate  missiles  were  fired  into  the  storm. 

The  radar  values  for  the  calibration  period  on  this  day  are  shown  in  the 
second  column  of  Table  5,  The  corresponding  Z\\  ^  values  (the  dll  equivalent  of 
Z\\'jj)  are  shown  in  the  third  column  and  were  obtained  from  one  of  the  two  dis- 
drometers  located  at  the  launch  site,  identified  as  "Disdrometer  E". 

Decibel -equivalent  values  of  the  system  calibration  constant,  that  is,  values 
of  jc  ,  were  determined  for  each  pair  of  the  Jj  and  Z\yn  values  as  follows: 

Equation  (85)  was  written  in  dll  form  as 

10  log  k  10  log «  +  10  log  I  -  10  log  Z...  ,  (94) 

n 

or,  in  the  dB  equivalent  form  defined  herein,  as 


o  + 


(95) 
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Table  5.  Example  of  Disdromster  Method  of  Radar  Calibration  for  17  February 
1972,  Disdrometer  B,  Decibel  Values,  Geometric  Averaging 


Time 

<dB) 

JL. 

<dB> 

fc-32 

<dB2) 

1458Z 

34.33 

73.0 

175.32 

+  0. 517 

0.27 

1502 

33.84 

70.8 

173. 61 

-  1. 193 

1.42 

1503 

32.42 

69.  3 

173. 53 

-  1.273 

1.62 

1504 

31. 14 

68.  5 

174.01 

-  0. 793 

0.  63 

1505 

30.36 

66.  3 

172. 59 

-  2.213 

4.90 

1 506 

28.61 

68.8 

176.84 

+  2.037 

4.  15 

1508 

27.57 

65.4 

174.48 

-  0.323 

0.  10 

1509 

31.00 

6f.  1 

172.75 

-  2.053 

4.21 

1510 

28.75 

66.  8 

174. 70 

-  0.  103 

0.01 

.511 

28.23 

67.  8 

176.22 

+  1.417 

2.01 

1512 

29.61 

67.0 

174.04 

-  0.763 

0.  58 

1513 

29.  97 

68.  1 

174.78 

-  0.023 

0.00 

151* 

30.81 

67.  1 

172.94 

-  1.863 

3.47 

1515 

29.  60 

66.  7 

173. 75 

-  1.053 

1.  11 

1516 

31.67 

68.8 

173.78 

-  1.023 

1.05 

1517 

29.65 

65.  3 

172.30 

-  2. 503 

6.27 

1518 

26.  67 

65.  3 

175.28 

»■  0.477 

0.23 

1519 

28.  11 

67.  3 

175.84 

+  1.037 

1.08 

1  520 

26.92 

67.  6 

177.33 

^  2.  527 

6.39 

1521 

28.94 

68.  3 

176.01 

r  1.207 

1.46 

1522 

28.  12 

69.  1 

177. 63 

*-  2.  827 

7.99 

1523 

27.77 

68.  9 

177. 78 

t  2.  977 

8.86 

1524 

28.  56 

69.  1 

177.  19 

+-  2.  387 

5.  70 

1 525 

32.98 

68.  3 

171. 97 

-  2. 833 

8.03 

1526 

31.04 

69.8 

175.41 

•  0. 607 

0.37 

S' 

n  25 

746.67 

1700.  5 

4370.08 

+  0.005 

71.91 

Mean 

29. 867 

68.  02 

174. 803 

'  0.000 

2.876 

Standard  Deviation  -  a  -  ±  1  996 

|K. 
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Then,  since  a  for  the  FPS- 18  radar,  see  Table  3  and  Eq.  (80),  has  the  value 
4.63  X  lO^mm^cm  m~^  giving^  136.05  dB,  Eq.  (95)  reduces  to 

k  -  136.  65  +  .1,  -  Z.  .  i90) 

^  L__J2j 

thereby  permitting  the  computations  cited  above. 

The  mean  value  of  ^  was  next  computed  for  the  n  samples  of  the  calibi’ation 
period.  This  mean  value, 


^  M 


_1_ 

IT 


i  n 


(97) 


was  assumed  to  be  the  value  of  the  system  calibration  constant  that  pertained  to  the 
particular  disdrometer  and  calibi’ation  period.  For  the  H  disdrometer  on 
17  February  1972,  the  Table  5  data  gave  the  value 


k  174,  80  dB 

i _ i  M 


(98) 


The  ^  values  for  the  other  disdrometers  and  days  of  the  1971-72  season 
were  computed  in  an  identical  manner.  The  resultant  values,  including  the  above, 
are  listed  in  Table  6. 

Additional  information  was  sought  in  the  1971-72  season  concerning  the 
standard  deviation  of  the  ^  ^  and  X\y ^  values,  also  regarding  the  correlation 

coefficient  (r)  and  the  standard  error  of  estimate  (S.  E. )  of  the  I  versus  Z;yn 
data.  The  values  of  these  parameters  are  also  shown  in  Table  6  for  the  different 
disdrometers  and  days  of  the  1971-72  season. 

The  only  parameters,  of  these  tabulated,  that  warrant  special  comment  are 
the  standard  deviation  of  ,k,  ,  that  is  a  ,  and  the  correlation  coefficient. 

•— 1  ’  Jij 

It  should  be  emphasized  that  the  a  values  of  Table  (>  are  not  indicative  of 

LIU 

the  uncertainties  of  radar  calibration  under  the  disdrometer  method.  The  values 
merely  indicate  the  uncertainty  or that  would  have  resulted  if  any  single  pair 
of  the  ^  and  Z\y^  values  had  been  used,  by  themselves,  for  radar  calibration 
purposes.  This,  of  course,  was  not  done  since  it  was  the  mean  value  ^ 
for  the  entire  sampling  period,  that  was  used  to  define  the  calibration  constant. 

The  correlation  coefficients  shown  in  Table  6  must  be  interpreted  with 
extreme  care.  The  correlation  coefficient  for  the  I  and  Z\yn  values  of  the  data 
sets  is  dependent  on  the  range  of  the  rainfall  rates  that  occurred  within  the 
calibration  periods.  Thus,  the  correlation  coefficients  will  (and  do)  have  larger 
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Table  G.  Values  of  the  System  Calibration  Constant,  Standard  Deviation  Values 
and  Correlation  Parameters  Obtained  by  the  Disdrometer  Method  for  the  FPS-18 
Radar  for  the  1971-72  Season 


Date  and 
Disdrometer 
Employed 

Number 
of  Paired 
Measure¬ 
ment 
Samples 

Calibration  Constant 
and  Standard  Deviation 

Correlation  Parameters, 
of  the  J, versus  Z...  values 
>— _J2i 

“M 

<dB) 

o 

,  K, 

<dB> 

°T 

ilj 

(dB) 

°Z 

A\\ 

L-  D. 
(dB) 

r 

*S.  E. 

(dll) 

12 

172.8 

1.8 

2.9 

2.  G 

0.791 

1.  6 

Disdrometer  B 

12 

173.4 

1.2 

2.9 

2.9 

0.  907 

1.2 

17  February  1972 

Disdrometer  A 

25 

175,8 

2.0 

1.7 

2.  1 

0.432 

1.9 

Disdrometer  B 

25 

174.8 

1.7 

1.7 

2.  1 

0.  590 

1.7 

17  March  1972 

Disdrometer  A 

45 

178.8 

2.3 

3 .  G 

3.9 

0.807 

2.3 

Disdrometer  B 

45 

1 79.  9 

2.9 

3.G 

4.7 

0.  78  G 

W3M 

22  March  1972 

Disdrometer  A 

28 

177.2 

4.  1 

3.8 

4.  G 

0,  553 

3.9 

Disdrometer  B 

28 

175. 1 

4.3 

3.8 

4.4 

0.453 

m a 

values  for  data  sets  in  which  there  is  a  substantial  range  of  the  rainrate  values 

(implying  an  also  substantial  range  of  the  separate  I  and  Z\v,x  values)  than  for 

' — *  i _ Lb 

data  sets  in  which  the  rainrates  are  more  or  less  constant  and  uniform  throughout 

the  calibration  periods.  To  state  this  another  way,  the  regression  line  for  the 

and  Z\vf,  points,  as  plotted  on  a  scatter  diagram,  will  tend  to  be  "flat”  (horizontal) 

for  a  uniform  rainrate  situation,  since  there  is  no  association  trend  of  the  Jj  and 

Z\\  n  values,  only  scatter.  *  The  correlation  coefficient  in  such  a  case  has  a  small 
i — —tii 

and  relatively  meaningless  value.  The  scatter  of  the  data  is  given  by  the  o  and 

i K  i 

S.  E.  values.  Conversely,  when  the  range  of  the  rainrates  during  the  calibration 
period  is  substantial,  the  regression  line  fitting  the  data  points  now  has  "slope" 
and  the  correlation  coefficient  lias  a  larger  value  than  before  (assuming  no  greater 
scatter  of  the  data).  Some  significance  can  be  attached  to  the  value  of  the  correla¬ 
tion  coefficient  but  interpretative  care  must  be  exercised,  because  the  value. 


-Strictly  speaking,  depending  on  which  parameter  is  correlated  with  the  other, 
the  regression  line  mav  be  horizontal,  vertical,  or  undefined;  hence,  cease  to 


in  part  and  degree,  is  still  dependent  on  the  particular  range  of  the  rainrate  values 

within  the  calibration  period.  The  scatter  of  the  data,  as  before,  is  given  by  the 

a  and  S.  E.  values. 
tL 

The  methods  used  to  obtain  the  system  calibration  constants  for  the  FPS-18 
radar  in  the  1971-72  season  have  been  explained  in  the  preceding  paragraphs.  It 
should  be  noted  that  the  constants  were  determined  by  averaging  decibel  values 
to  obtain  geometric  means. 

A  question  may  be  raised  concerning  the  differences  in  the  means  and 
standard  deviations  that  would  have  resulted  if  the  actual,  or  absolute,  values  of 
I  and  had  been  averaged  to  obtain  arithmetic  means.  These  differences 

were  explored  for  the  data  of  the  II  disdrometer  for  17  February  1972,  the  same 
data  set  used  previously  for  illustration.  The  results  are  shown  in  Table  7. 

A  comparison  of  Tables  5  and  7  shows  that  the  value  determined  by 
arithmetic  averaging  has  the  dB-equivalent  value  of  173.  14,  which  differs  by 
^0.  34  dB  from  the  value  of  T  V  5.  It  is  also  seen  that  the  standard  devia¬ 
tion  value  of  Table  7  in  terms  of  its  dB-equivalent,  computed  as  indicated  at  the 
bottom  of  the  table,  is  ±1.  835  dll  as  opposed  to  a  value,  from  Table  5,  of 
tl.  090  dB,  a  difference  of  ±0.  139  dll. 

This  indicates  the  magnitude  of  the  differences  between  the  values  of  the 
calibration  constants  that  might  be  anticipated  from  arithmetic  averaging  as 
opposed  to  geometric  averaging.  The  differences  are  relatively  small  but  are 
nevertheless  significant,  particularly  in  view  of  the  SAMS  objective  of  minimizing 
all  possible  sources  of  error  and  uncertainty.  The  question  of  which  method  of 
averaging  is  preferred  is  really  an  open  question  since  many  radar  parameters, 
such  as  J_,  ,  are  measured  as  dll  averages;  complex  problems  exist  concerning 
the  relationship  of  these  averages,  and  standard  deviations,  to  the  absolute 
values.  For  example,  see  Lhermitte,'*0  Lhermitte  and  Kessler,  31  and  Austin  and 
Schaffner.  18 

To  summarize  the  discussions  of  the  present  section  concerning  the  dis¬ 
drometer  method  of  calibration,  reference  is  made  to  Table  8.  The  system 
calibration  constants  for  the  FPS-18  radar  in  the  1971-72  season  are  listed  in  this 
table  in  terms  of  their  absolute  values,  rather  than  the  dB  values  of  Table  8. 


30.  Lhermitte,  It,  M.  (1903)  Motions  of  Scatterers  and  the  Variance  of  the 
Mean  Intensity  of  Weather  Radar  Signals.  Research  Report  03-57,  Sperry  Rand 
Research  Center,  Sudbury,  Massachusetts. 

31.  Lhermitte,  R.M.,  and  Kessler,  E.  (1906)  Estimation  of  the  Average 
Intensity  of  Precipitation  Targets.  Proceedings  12th  Conference  on  Radar 
Meteorology,  17-20  October  1906,  Norman,  Oklahoma. 
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Table  7.  Example  of  Disdrometer  Method  of  Radar  Calibration  for  17  February 
1972,  Disdrometer  B,  Absolute  Values,  Arithmetic  Averaging 


Standard  Deviation  -  a  ±  1,305  <  10  *7 

K 

i T  +  a  4.  572  x  10^ 7 

K 

176.60  dR 

7  -  a  -  1. 962  <  1017 

K 

172,  93  dB 

Difference 

3.  67  dB  i-2'  odH'  t  1. 83  5  dB 

GO 

Table  8.  Absolute  Values  of  the  System  Calibration  Constant  k  and  Associated 
Radar  Parameters  Determined  by  the  Disdrometer  Method  for  the  FPS-18  Radar 
for  the  1970-72  Season 


Corresponding  Values  of 

Date  and 

Calibration 

1. 

k 

C 

1 

Disdrometer 

Constant 

4  -•> 

employed 

K 

m\\(.\M)  cm  " 

(NM)“t 

m  m\V 

*15  March  1971 

1 26  March  1971 

l  n known 

I'nknown 

0.0533 

1.259  s  10*4 

5.14  • 

10:> 

3  February  1972 

Disdrometer  A 

1  7 

1.91  «'  10‘ 

0. 0296 

0.699  v  10'4 

3.85  - 

105 

Disdrometer  11 

17 

2. 19  '  10 

0.0340 

0.802  •  io‘4 

4.42  ■ 

io5 

17  February  1972 

Disdrometer  A 

3.80  <  1017 

0. 0589 

1.390  >  10-4 

7.  66  \ 

io5 

Disdromrler  11 

1 7 

3.02  \  10l 1 

0.0468 

1.  104  <  10‘4 

6.08  •' 

10:> 

17  March  1972 

Disdrometer  A 

17 

7. 50  v  10A  * 

0.  1177 

2.778  v  10'4 

15.30  • 

10'’ 

Disdrometer  H 

1  7 

9.77  v  101 

0. 1515 

3.575  ■  10"4 

1  9,  70  s 

ior’ 

22  March  1972 

Disdrometer  A 

17 

6.76  v  10x 

0.  1048 

2.473  \  10'4 

13,62  ' 

10-’ 

Disdrometer  R 

1  7 

3.71  -■  10 

0. 057  5 

-4 

1.357  ■  10  * 

7.48  v 

10  5 

Tiie  k  values  for  the  1970-71  seasons  arc  unknown  because  a  variable 
coupling  attenuator  was  used  with  the  video  integrator.  The  other  values 
listed,  however,  are  values  determined  from  disdrometer  calibration. 


Also  listed  are  the  values  of  the  system  loss  factor, 


see  Kqs.  ,7b)  and  (80),  the  sphere  calibration  constant, 


(99) 


(100) 


see  Kc|.  (78),  and  the  calibrated  constant  C  of  the  radar  equation  for  hydrometeor 
targets,  computed  from  Kq. (83).  The  C’  and  c  values  which  enter  Kq.  (99)  are 
given  in  Table  4. 
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The  k.  L,  k  and  C  values  for  the  1970-71  season  have  also  been  included  at 
'  *  s 

the  top  of  Table  8. 


The  C.  til  i  brut  ion  Constant*  of  \pplieution  und  Kslimute*  of  the  Possible 
I  neerluinlies  of  Calibration 

Since  two  different  methods  of  radar  calibration  were  employed  for  the 
FPS-18  radar  in  the  1970-71  and  1971-72  seasons,  there  was  the  obvious  question 
of  which  of  the  calibration-constant  values  should  actually  be  applied  to  the  radar 
data  for  the  SAMS  missile  trajectories.  This  question  will  be  considered  in  the 
present  section  and  we  will  begin  by  discussing  the  relative  merits  and  dis¬ 
advantages  of  the  two  calibration  methods.  Comments  will  also  be  made  about  the 
possible  uncertainties  of  calibration. 

The  sphere  method  of  calibration  lias  the  rather  obvious  advantage  of  simplicity. 
The  value  of  the  calibration  constant  is  readily  obtained  from  the  spherical 
targets  tracked  by  the  radar.  The  method  lias  one  severe  disadvantage,  however, 
in  that  the  calibration  is  almost  invariably  performed  on  a  day  other  than  that  of 
missile-launch  operations.  Thus,  if  the  calibration  constant  obtained  from  sphere 
calibration  is  applied  to  the  radar  data  along  the  missile  trajectories,  it  is 
implicitly  assumed  that  the  radar  on  the  launch  day  was  operating  the  same  as  it 
was  on  the  calibration  day.  This  may  or  may  not  be  true. 

There  arc  several  reasons  why  sphere  calibrations  are  performed  on  days 
other  than  the  launch  days.  First  of  all,  the  calibration  operations  require  about  a 
one  to  two  hour  period  of  radar  time,  together  with  the  time  of  the  personnel 
operating  the  radar.  On  a  day  of  launch  operations,  since  the  radar  is  used  for 
launch  planning  purposes,  it  is  difficult  to  schedule  the  time  needed  for  sphere 
calibration  and  also  tolerate  the  diversion  of  personnel  to  the  task.  More 
fundamental  though,  it  has  been  found  that  if  a  sphere  calibration  is  attempted  on 
the  day  of  launch  operations— which  is  usually  a  stormy  day,  almost  by  definition  — 
the  balloon  or  parachute  suspending  the  sphere  becomes  wet  from  the  precipitation 
and  gives  echo  return  on  the  radar  that  destroys  the  validity  of  the  calibration 
measurements. 

It  is  difficult  to  estimate  the  uncertainties  of  radar  calibration  associated  with 
the  sphere  method.  To  enable  such  estimates,  comparison  calibrations  would 
have  to  be  performed  on  one  si.  gle  day.  Such  comparisons,  to  the  author's  know¬ 
ledge,  have  not  been  made. 

The  disdrometer  method  of  radar  calibration,  in  contrast  to  the  sphere, 
method,  has  the  advantage  of  being  conducted  on  the  same  day  as  the  launch  opera¬ 
tions,  in  the  same  storm  conditions  (or  very  similar)  encountered  by  the  missiles. 
The  prime  disadvantage,  at  least  in  the  past,  is  that  the  radar  and  disdrometer 
measurements  required  for  calibration,  also  the  computations  and  analyses,  are 


considerable  and  time-consuming  such  that  the  values  of  the  calibration  constants 
could  not  be  obtained  until  some  weeks  or  months  following  the  launch  day. 

There  are  several  other  problems  and  questions  concerning  the  disdrometer 
method  of  calibration  that  are  worthy  of  mention.  There  is  the  problem  of  the 
distance  separation  of  the  places  of  disdrometer  and  radar  measurement.  This 
separation  undoubtedly  contributes  importantly  to  the  rather  large  scatter  of 
the  comparative  data.  The  calibration  constant  itself,  which  is  a  mean  obtained 
for  the  entire  period  of  common  sampling,  should  in  theory  be  reasonably  valid, 
providing  that  the  spatial  inhomogeneities  of  rainfall  are  more-or-less  randomly 
distributed  throughout  the  Wallops  storms.  But  whether  this  is  actually  the  case 
and,  if  not,  whether  "separation  bias"  exists  within  the  data  are  unanswered 
questions. 

There  are  also  questions  concerning  the  proper  sensitivity  tests,  calibration- 
procedures,  and  siting  and  exposure  methods  that  should  be  used  with  the  dis¬ 
drometer  instruments  to  insure  their  accurate,  optimum  performance.  In  fact, 
quantitative  information  is  lacking  as  to  what  really  constitutes  optimum  perfor¬ 
mance  under  the  variety  of  storm  conditions  encountered  at  Wallops  Island. 
Problems  likewise  exist  concerning  the  large  disparity  between  the  sample  volumes 
of  the  disdrometers  (of  the  order  of  0.3  to  3m'*)  and  those  of  the  radar  (of  the  order 
of  105  to  107  m*). 

In  view  of  these  several  questions  and  problems,  it  is  as  difficult  to  assess  the 
possible  uncertainties  of  calibration  for  the  disdrometer  method,  as  for  the  sphere 
method.  It  might  be  argued  that  the  uncertainties  for  the  disdrometer  method 
might  be  given  by  the  differences  in  the  calibration  constants  computed  from  the 
data  of  two  adjacent  disdrometers  on  the  same  calibration  day.  This  might  be  true, 
to  a  first  approximation,  but  it  should  also  be  pointed  out  that  the  data  sets  for  the 
two  disdrometers  are  not  truly  independent,  since  the  same  radar  signal  values 
were  employed  commonly  in  both. 

Because  of  the  difficulties  with  the  two  calibration  methods  indicated  above, 
it  was  decided,  at  least  for  the  1971-72  season,  that  the  final  calibration  values 
to  be  applied  to  the  radar  measurement  data  for  the  missile  trajectories  would  be 
the  simple  average  of  (1)  the  sphere  value,  (2)  the  value  for  disdrometer  A  and 
(3)  the  value  for  disdrometer  B.  This  was  an  arbitrary  decision  that  effectively 
weighted  the  values  of  the  disdiometer  method  twice  those  of  the  sphere  method. 

These  final  calibration  constants,  thus  computed  for  the  FPS-18  radar  for 
the  1971-72  seasor,  are  shown  in  Table  9.  Also  shown  are  the  values  of  the 
calibration  constants  that  were  applied  to  the  FPS-18,  the  X  band  and  the  Spandar 
radars  during  the  other  SAMS  seasons. 

-This  will  be  ameliorated  in  the  future  by  greater  automation  of  the  dis¬ 
drometer  equipment  and  recording,  plus  a  more  extensive  use  of  computers. 
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Table  9.  Final  Values  of  the  System  Calibration  Constant  and  Associated  Parameters  Used 
for  the  SAMS  Radar  Measurements  of  the  Different  Seasons 
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•'•Unknown  because  variable  coupling  attenuator  used. 

-See  Table  4  footnote  (f)  for  explanation  of  the  two  coupling  relations  used  with  the  Spandar  radar  in  the 
1972-73  season. 


As  a  last  comment  concerning  the  uncertainties  of  calibration,  it  is  suggested 
that  these,  for  the  FPS-18  radar  in  the  1971-72  season,  might  be  indicated  to  a 
first,  gross  approximation  by  the  differences  in  the  two  disdrometer  values  and 
the  sphere  value  that  were  used  to  compute  the  Table  9  averages.  The  difference 
spread  of  the  values  on  the  four  days  of  the  1971-72  season,  see  Tables  4  and  8, 
ranged  from  0.  6  to  4.  6  dB  (equivalent)  and  it  may  be  stated  that  the  standard  error 
of  estimate  varied  from  ±0.4  dB  to  ±3.3  dB.  The  calibration  uncertainties  are 
probably  of  this  order. 


3.0  Review  of  the  Trujerlorv  Measurement  Procedure* 

The  discursive  continuity  of  the  report  has  been  interrupted,  in  the  last  three 
sections,  by  the  necessity  of  describing  the  radar  calibration  methods  in  some 
detail.  Therefore,  it  would  seem  beneficial,  at  this  point,  to  review  and  con¬ 
solidate  the  material  presented  thus  far  and  to  summarize  it  in  the  context  of  its 
direct  application  to  the  SAMS  measurement  objectives. 

In  review,  then,  we  might  first  reiterate  that  the  SAMS  missiles  were  fired 
in  the  14G°  azimuth  direction  from  launch  pad  zero  at  Wallops  Island  and  that  the 
nominal,  average  elevation  angle  of  the  trajectories  from  the  launch  pad  was  about 
26°  to  31°. 

Radar  data  for  the  missile  trajectories  were  obtained  from  the  RAUF  radars 
which  were  located  1.8  NM  to  the  NW  of  the  launch  pad.  The  data  were  obtained 
primarily  in  RBI  scan  in  the  140°  azimuth  direction  from  the  radars,  which 
very  nearly  corresponded  within  1/2°  or  so,  to  the  140°  direction  from  the  launch 
pad.  Video  integration  was  used  and  the  I  signal  values  for  each  of  the  RIII  scans 
were  photographically  recorded,  in  the  1970-71  and  1971-72  seasons  (as  gray  scale 
regions  of  shading  contrast),  and  were  recorded  on  video  tape  for  computer  proces¬ 
sing,  in  the  1972-73  season. 

Once  the  missile  trajectories  had  been  established  from  the  NASA,  postflight 
tracking-data,  the  I  values  along  the  course  of  the  trajectories  were  determined 
from  the  radar  RHI  photographs,  or  computer  products,  that  pertained  to  the  launch 
times  of  the  missiles.  These  I  values  then  became  the  initial,  input  data  for  the 
computation  of  the  desired  end-product  values  of  the  reflectivity  factors  Z^  for 
water  hydrometeors,  and  Zj  for  ice  hydrometeors. 

To  obtain  quantitative  values  of  Z^.  and  Z^,  two  things  had  to  be  accomplished. 
Both  were  essential  and  the  order  of  mention  is  unimportant.  First,  the  radar  had 
to  be  calibrated.  This  was  accomplished  by  the  two  methods  described  previously. 
The  particular  calibration  constant  of  computational  interest  was  the  system 
calibration  constant,  k,  the  final  values  of  which  arc  shown  in  Table  9.  Second, 
independent  knowledge  was  required  for  the  Wallops  storms  concerning  the  precise 
altitude  bounds  of  the  melting  zone,  which  separated  the  region  of  ice  hvdrometeors 
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aloft  (the  ice  crystal  and  snow  particles)  from  the  rain  region  below.  Knowledge 
of  the  storm  top  altitude  was  also  required,  since  the  RARF  radars,  except 
under  particular  conditions  (of  near  vertical  incidence)  are  insufficiently  sensitive 
to  detect  the  storm  top  boundary.  This  knowledge  was  obtained  from  aircraft 
measurements  and  observations,  as  will  be  elaborated  in  a  subsequent  report. 

With  the  radars  calibrated,  the  water  and  ice  regions  defined,  and  the  storm 
top  altitude  established,  the  Z^-  values  for  the  water  region  and  particular  radar 
were  computed  from  the  T  data  using  Eq.  (72),  that  is. 


1.08  X4  I  v  .  n12  mmC 
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and  the  Zj  values  for  the  ice  region  and  particular  radar  were  computed  from 
Eq.  (73),  that  is, 


^  4. 78  A4  I  v  ,„12  mmG 

ZI  5 -  X  10 !T~ 

r r  k  m 


For  example,  for  the  FPS-18  radar  on  17  February  1972— a  radar  and  a  day 
previously  used  in  illustration  — Eqs.  (721  and  (73)  immediately  reduce  to  (see 
Tables  3  and  9) 
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Withtn  the  melting  zone,  as  noted  in  Section  3.  l,  the  values  of  the  reflectivity 
factor  are  indeterminate.  The  values  of  both  'V  and  Zj  were  computed  for  the 
segments  of  the  missile  trajectories  that  passed  through  this  zone.  The  Z^. 
values  are  indicative  of  the  small-value  limits  of  the  reflectivity  factor  across  the 
zone  and  the  Zj  values  provide  some  indication  of  the  reflectivity  factor  values 
which  are  equaled  or  exceeded  within  the  zone.  The  explanation  of  this  last 
statement  is  that  snowflakes,  which  fall  into  the  melting  zone  and  acquire  a  thin 
water  coating  without  appreciable  size  change,  can  produce  radar  returns  greatly 
exceeding  those  of  dry  snowflakes  of  equal  size. 


Another  special-case  situation  was  also  mentioned  previously.  This  was  the 
case  of  supercooled  clouds,  which  exist  in  water  form  at  temperatures  well  below 
freezing.  It  may  be  stated  that  no  analytical  problems  existed  with  the  radar 
data  regarding  such  clouds,  or  regarding  water  clouds  in  general,  either  warm  or 
supercooled,  since  the  RARF  radars  are  presently  incapable  of  detecting  water 
droplets  in  the  cloud-size  range,  SAMS  information  about  water  clouds  had  to  be 
acquired  necessarily  from  aircraft  measurements. 


1.  RKSCRII’TION  \M>  PHORHKSS  HKI'OHT  OK  IRK  SW1S  IUDVR  WK  VSI  RKMKNTS  OK  TIIK 

IMI  KKRKVr  SKVSONS 

The  general  nature  of  the  SAMS  operations  and  radar  measurements  of  the 
different  seasons  to  date  will  be  described  in  this  section  and  comments  will  be 
made  about  plans  for  the  1973-74  season. 

These  descriptions  will  serve  two  purposes.  They  will  indicate  the  particular 
types  of  radar  measurements  that  were  made  in  the  different  years  and  they  will 
provide  a  progress  report  to  show  how  the  specificity  and  accuracy  of  measure¬ 
ment  have  been  improved,  season  by  season,  to  better  satisfy  the  SAMS  overall 
objectives. 

It  should  be  noted  that  missiles  of  the  Terrier-Recruit  type  were  used 
throughout  all  of  the  previous  seasons.  These  missiles  and  flights  have  been 
reported  by  Cole  and  Robinett,  and  Cole,  Church  and  Marshall. 

1. 1  The  1170-7 1  Season 

The  SAMS  program  at  Wallops  Island  was  in  its  infancy  in  the  1970-71 
season.  AFCRL  had  no  prior,  extensive  experience  with  radar  measurements  of 
the  type  required  for  the  program.  The  Wallops  Island  radars  and  video  integrator 
were  untested  and  the  calibration,  measurement  and  analytical  procedures  had  to 
be  developed  as  the  season  progressed. 

Both  the  FPS-18  and  X  Band  radars  were  used  to  obtain  data  along  the  missile 
trajectories.  The  FPS-18  data  consisted  of  RHI  and  A  scope  photographs  of  the 
RIfl  sweeps  made  in  the  launch  direction.  The  X  band  data  consisted  of  graphical 
records  of  received  power  versus  range  that  were  acquired  for  a  fixed  elevation 
angle  corresponding  approximately  to  the  mid -trajectory  angle  of  the  missile,  as 
measured  from  the  origin  at  the  radar  site.  The  records  were  nicknamed  the 
"Bloody  Reds'  ,  since  the  traces  were  highly  variable  and  copious  amounts  of  red 
ink  were  expended  in  the  recording  process. 


The  X  band  data  for  this  season  were  analyzed  by  the  Applied  Physics 

Laboratory  of  Johns  Hopkins  University.  The  results  have  been  reported  by 
a  o 

Meyer. 

The  calibi'ation  constants  for  the  FPS-18  radar  were  obtained  by  both  the 
sphere  and  disdronieter  methods.  Those  for  the  X  band  radar  were  obtained  by 
the  sphere  method.  The  values  of  the  constants,  as  best  known  to  the  author,  are 
shown  in  Tables  4  and  8. 

Video  integration  was  used  just  with  the  FPS-18  radar.  The  integrator  pro¬ 
vided  a  running-mean  average  of  the  signal  ( (  1 1 )  values  which  was  "weighted"  in 
an  exponential  fashion.  The  averaging  was  performed  for  single  range  cells. 

The  integrator  also  had  an  "attenuator"  which,  in  effect,  permitted  the  values  of 
the  coupling  constant,  see  Eq.  (81),  to  be  changed  at  will. 

The  averaging  factor  y  for  this  integrator  and  season,  sec  Eq.  (44),  was 

assumed  to  have  the  value  0.398  (or  -4.0  dB),  a  value  obtained  from  the  work  of 
1 8 

Austin  and  Schaffner. 

It  is  the  author's  understanding  that  the  video  integrator  of  the  1970-71  season 
was  deficient  in  two  primary  respects.  First,  it  failed  to  provide  integration  for 
a  sufficient  number  of  independent  samples,  which,  combined  with  the  exponential 
weighting  of  the  averages,  resulted  in  the  I  signal  output  being  somewhat 
"jittery”.  This  tended  to  degrade  the  quality  of  the  gray-scale  photographs 
showing  the  signal  levels  within  the  Wallops  storms.  The  integrator  also  failed  to 
accomplish  the  "range  normalization  function",  see  Section  2,7.2,  with  complete 
fidelity,  such  that  range  corrections  had  to  be  applied  in  certain  instances. 

The  data  for  each  missile  were  obtained  from  the  Hill  photographs  taken 
at  the  launch  times.  The  course  of  the  given  missile  trajectory  across  tlu  photo¬ 
graph  was  known  and  the  Jj  signal  values  could  be  accurately  determined  for 
each  p  >int  where  the  trajectory  intersected  with  a  boundary  of  the  gray-scale 
regions  shown  on  the  photographs.  The  Jj  values  between  points,  that  is  the 
values  along  the  trajectories  within  the  photographic  gray-regions  themselves, 
could  only  be  inferred,  or  estimated,  by  interpolation. 

The  trajectory  values  of  Z^.  and  Zj  were  determined  from  the  Jj  data  follow  ¬ 
ing  the  procedures  outlined  in  Section  3.8.  Information  about  the  altitudes  of  the 
malting  zone  and  storm  top  level  was  acquired  by  the  AFCHI.,  C-130,  cloud- 
physics  aircraft. 


32.  Meyer,  J.  II.  (1971)  Precipitation  Hates  and  Liquid  W  ater  Content 
Measured  by  the  JAI’NA  X-hand  Hadar  for  the  ltain  Erosion  Test.  29  March  1971, 
Informal  Report,  Applied  Phvsics  Laboratory,  Johns  Hopkins  University, 
BPD714-1 9,  24  June  1971. 


G8 


1.2  The  1 ‘>71-72  Season 


The  FPS-18  radar  was  used  as  the  prii  ~y  radar,  as  in  the  previous  season, 
and  RHI  and  A  scope  data  were  obtained  for  e  missile  trajectories  in  a  manner 
similar  to  the  first  season,  but  with  certain  improvements.  The  Spandar  radar, 
rather  than  the  X  band,  was  used  to  secure  the  "Bloody  Red"  records  of  Pr  versus 
range,  for  the  trajectories. 

The  FPS-18  radar  was  calibrated  by  both  the  sphere  and  disdrometer  methods, 
as  described  in  Sections  3.3,  3.4  and  3.  5  herein.  The  Spandar  radar  was  cali¬ 
brated  by  the  sphere  method  only.  The  calibration  constants  of  final  application 
are  listed  in  Table  9. 

A  video  integrator  of  new  design  was  utilized  in  the  1971-72  season.  This 
provided  logarithmic,  running-mean  averaging  of  the  instantaneous  signal  values 
and  its  versatility,  in  terms  of  the  selectable  choices  of  integration  period,  i. umber 
of  range  cells  included,  number  of  signal  level  outputs,  and  output  format,  was  a 
considerable  improvement  over  the  integrator  used  in  the  previous  season.  The 
averaging  factor  ->  for  this  integrator  and  season,  was  assumed  to  be  0.  a37 

1 B 

(or  -2.  7  dB),  a  value  again  obtained  from  the  work  of  Austin  and  Schaffner. 

The  quality  of  the  grayscale,  RHI  photographs  acquired  from  the  FPS-18 
radar  was  enhanced  by  the  employment  of  the  new  integrator.  A  larger  number  of 
shading  differentiations  of  signal  level  was  provided  on  the  photographs,  than 
formerly,  and  the  contrast  definition  of  the  gray-shade  regions  provided  better 
and  more  accurate  picture  resolution.  The  photography  was  accomplished  using 
Polaroid  film,  which  was  developed  for  viewing  some  30  seconds  follow  ing  the 
termination  of  the  RHI  sweep.  Thus,  the  photographs  could  be  used  for  real-time 
planning  purposes  and  were  also  available  for  later  analyses. 

The  values  for  the  missile  trajectories  were  obtained  from  the  launchtime, 
RHI  photographs  in  the  same  manner  as  in  the  previous  season.  The  values  were 
determined  for  each  point  along  the  trajectory  which  intersected  the  boundaries 
of  the  gray-scale  regions  shown  on  the  photographs.  The  trajectory  values  between 
points  were  inferred  by  interpolation. 

The  Z^,  and  Zj  values  for  the  trajectories  were  established  from  the  data 
using  Eqs.  (72)  and  (73),  which  were  programmed  for  computer  solution  and 
plotting.  Information  about  the  altitudes  of  the  boundaries  of  the  melting  zone  and 
storm  top  was  secured  by  the  AFCRL,  C-130  aircraft  and/or  by  an  Aztec  aircraft 
flown  contractually  by  Meteorology  Research,  Inc.  (MRI). 

The  trajectory  data  acquired  from  the  Spandar  radar  were  never  fully- 
analyzed  for  the  1971-72  season.  The  analyses,  begun  by  APL,  were  time- 
consuming  and  there  was  an  apparent  general  realization  among  the  AFCRI.  and 
API.  personnel,  which  was  never  verbalized,  that  the  Spandar  data  were  not 
likely  to  provide  much  of  additional  value  to  the  information  already  obtained 
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from  the  FPS-18  photographs.  The  analyses  were  discontinued,  essentially  by 
default  of  interest. 

1.3  The  1172*73  Season 

The  FPS-18  radar  was  also  used  as  the  primary  radar  in  the  1972-73 
season.  The  "Bloody  Red"  recordings  %vere  discontinued  and  the  Spandar  radar 
was  diverted  to  operation  in  the  PPI  mode,  to  secure  PPI  and  sector-scan 
photographs  of  polaroid  type. 

Video  integrators  were  used  with  both  of  these  radars  which  were  of  the 
same  type  as  the  previous  season.  The  integrator  for  the  FPS-18  radar  pro¬ 
vided  contoured  Kill  information,  as  before,  and  the  Spandar  radar  provided 
contoured  PPI  and  sector-scan  data. 

A  major  change  was  effected  in  the  1972-73  season  concerning  the  methods 
of  recording  the  primary  data  from  the  FPS-18  radar.  Although  photographic 
recording  wat  still  used  for  real-time  planning  purposes,  the  integrator  output 
for  the  RHI  sweeps  was  also  recorded  on  video  tape.  These  tapes  were  processed 
by  computer  using  programs  designed  by  the  Weather  Radar  Branch  of  AFCRL. 
Various  kinds  of  information,  in  RHI  format,  were  provided  and  greatly  simpli¬ 
fied  the  analytic  procedures  of  determining  the  signal  values  for  the  missile 
trajectories.  The  accuracy  and  resolution  of  the  analyses  were  also  enhanced, 
relative  to  the  photographic  procedures  used  previously. 

The  method  of  determining  the  Jj  values  along  the  missile  trajectories  is 
illustrated  in  Figure  2,  A  rectangular  gridwork  is  shown  which  encompasses  a 
part  of  the  range -altitude  section  of  one  of  the  RHI  sweeps  of  the  FPS-18 
antenna.  (The  date  and  time  are  unimportant,  since  the  illustration  is  schematic 
only. )  It  is  seen  that  a  typed  letter,  number,  or  other  symbol  appears  within 
each  of  the  boxes.  These  provide  coded,  computer-output  information  about  the 
average  ^  signal  values  that  existed  within  the  range -altitude  boundaries  of  the 
boxes.  The  box  dimensions  are  f>00  ft  in  horizontal  extent  by  180  ft  in  vertical 
extent,  dimensions  that  are  dictated  in  part  by  the  letter  and  line  spacing  of  the 
teletype  printer  of  the  computer. 

The  beginning,  lower  portion  of  a  simulated  missile  trajectory  has  been 
drawn  on  the  diagram  which  cuts  through  and  between  the  various  boxes  along  its 
course.  It  is  seen  that,  with  knowledge  of  the  symbol  code,  the  Jj  %'alues  along 
the  trajectory  can  be  readily  established  as  a  function  of  distance.  Some  inter¬ 
polation  is  required  but  this  poses  no  particular  problem. 

This  procedure  was  used  to  secure  the  trajectory  data  for  the  1972-73 
season.  It  will  also  be  used  in  the  1973-74  season. 
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Radar  calibration  in  the  1072-73  season  was  accomplished  only  by  the  sphere 
method.  Data  were  acquired  to  permit  calibration  by  the  disdrometer  method  but, 
as  mentioned  earlier,  preliminary  analyses  revealed  discrepancies  between  the 
sphere  and  disdrometer  values  for  the  day  of  27  February  1073,  which  have  not 
yet  been  resolved. 

Determination  of  the  and  Zj  values  from  the  I  trajectory  data  was  ac¬ 
complished  as  described  in  Section  3.(>,  The  altitudes  of  the  melting  zone  boundaries 
and  storm  top  were  obtained  from  the  AFCR1.,  C-130  aircraft  and  or  from  a 
Navajo  aircraft  flown  contractually  by  MRI. 

1. 1  Nan*  for  the  I1!?;!-?  I  Sruson 

New,  three-stage,  3' AT  lilt  (Talos-Tcrrier-llecruit)  missiles  will  be  flown 
in  the  1973-74  season.  The  trajectories  of  these  missiles  will  be  lower  (more 
horizontal)  than  those  of  the  previous  two -stage  missiles  and  the  total  path  distance 
through  the  Wallops  storms  will  be  extended  appreciably.  This  imposes  the  require¬ 
ment  that  the  radar  measurements  for  the  new  trajectories  will  have  to  be  extended 
to  greater  radar  ranges,  than  for  the  old  trajectories.  It  also  imposes  a  radar 
detectability  problem,  since  the  detection  capabilities  of  the  FPS-18  radar  for  the 
small-snow  particles  and  ice  crystals  in  the  uppermost  portions  of  the  storms  is 
marginal,  at  best,  even  with  the  prior  trajectories. 

Computations  indicate  that  a  5  dll  increase  in  the  sensitivity  of  the  FPS-18 
radar  is  needed  to  obtain  measurement  data  for  the  new  trajectories  that  are 
equivalent,  in  terms  of  minimum  detectable  values  of  the  reflectivity  factor,  to 
the  trajectory  data  obtained  previously.  A  10  dll  increase  is  needed  to  ensure  the 
consistent  detection  of  the  small -snow  particles  and  ice  crystals  at  the  storm  top 
extremities  of  the  new  trajectories. 

To  meet  these  needs,  the  Weather  Radar  R ranch  of  AFCRI.  plans  to  increase 
the  sensitivity  of  the  FPS-18  radar  by  about  12  dll.  Three  dll  will  be  obtained  by 
doubling  the  pulse  length  of  the  radar.  The  other  9  dll  improvement  will  be 
achieved  by  compensating  for  the  non-logarithmic  performance  of  the  radar  receiver 
at  power  levels  below  about  -90  dRm.  Such  compensation,  made  electronically  or 
by  computer  processing  of  data,  will  permit  the  integrated  received  power  P  , 
see  Eqs.  (42)  and  (30),  to  be  accurately  determined  down  to  levels  of  approximately 
-101  dRm  rather  than  -92.3  dRm,  as  at  present,  in  the  1972-73  season. 

The  Spandar  radar  in  the  1973-74  season  will  be  used  for  three  purposes.  It 
will  be  operated  to  secure  contoured  PPI  photographs  of  the  W  allops  storms  at 
about  five  minute  intervals  during  periods  before  and  after  the  launch  time  of  the 
missile.  Retween  the  times  of  these  photographs,  the  antenna  will  be  |K>inted  over 
the  disdrometer  sites  to  obtain  calibration  data.  At  the  launch  time  of  (lie  missile, 
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the  radar  will  be  operated  in  the  R III  mode  and  five  Kill  photographs,  plus 
recorded  video  data,  will  be  obtained  for  a  5°  azimuth  sector  symmetrical  about 
the  azimuth  angle  of  the  missile  trajectory. 

The  Spandar  radar,  for  some  or  all  of  these  measurements,  will  be  operated 
in  a  frequency-diversity  mode  in  which  the  wavelength  of  transmission  and 
reception  is  changed  from  pulse  to  pulse  of  the  radar,  sequentially  increased  then 
decreased,  by  an  amount  slightly  execding  the  bandwidth  of  the  radar  receiver. 

This  mode  of  operation  is  advantageous  in  that  there  is  then  no  need,  in  the  video 
integration,  to  "wait  for"  the  hydrometeors  to  redistribute  themselves  naturally 
(through  turbulent  and  differential  fall -velocity  mechanisms)  to  secure  independent 
samples.  Rather,  independent,  non-autoeorrelated  samples  of  the  hydrometeors 
are  obtained  for  each  pulse  of  the  radar.  Hence,  with  frequency  diversity  the 
Atlas  equation,  Eq.  (41)  specifying  the  "reshuffling  time"  of  the  hydrometeors  no 
longer  pertains,  and  the  number  of  independent  samples  per  integration  period  for 
the  case  of  a  stationary  antenna  is  given  by 

n.  Atj  prf  ,  (103) 

instead  of  by  Eq.  (A7)  of  Appendix  A,  as  formerly  provided.  Likewise,  the  other 
equations  of  Appendix  A,  which  contain  n.  as  an  independent  variable,  namely 
Eqs.  (A13),  (A14),  (A16),  (A18),  (A19),  (A24)  and  (A26),  are  all  modified 
accordingly. 

The  representativeness  of  the  Spandar  data  acquired  with  frequency  diversity 
will  be  considerably  improved,  as  compared  to  previous  data.  The  number  of 
independent  samples  per  integration  period  will  be  increased  by  the  approximate 
factor  tx  prf.  Consequently,  some  17  times  more  independent  samples  will  be 
included  in  the  integration  than  formerly  (for  a  prf  of  960  pps,  which  is  commonly 
used). 

The  Jj  data  from  the  EPS- 18  radar  for  the  missile  trajectories  will  be  analyzed 
in  a  manner  similar  to  that  described  for  the  1972-73  season.  The  trajectory  data 
from  the  Spandar  radar  will  be  analyzed  contractually  by  API.. 

The  altitudes  of  the  melting  zone  boundaries  and  storm  top  will  be  determined 
by  the  AFCRL,  C-130A  aircraft  and/or  by  Navajo  or  Citation  aircraft  operated  by 
MRI,  A  second  AFCRL,  C-130E  aircraft  should  be  available  for  the  SAMS 
observation-measurement  program  about  March  1974, 
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Appendix  A 

Other  Considerations  Involving  Pulse  and  Video  Integration 


I.  INTItOIN  CTION 

When  radar  data  are  acquired  by  use  of  a  pulse  or  video  integrator,  it  is 
essential,  as  noted  in  Section  2. 4,  to  have  quantitative  information  concerning 
(1)  the  temporal  and  spatial  resolution  of  the  measurements,  (2)  the  independence 
and  representativeness  of  the  hydrometeor  samples  included  in  the  integration,  and 
(3)  the  variance  of  the  measurement  values  about  the  probable  "true  values". 

The  above  topics  are  discussed  herein.  The  discussion  provides  background 
information  regarding  radar  measurements  obtained  with  pulse  or  video  integrators 
which  is  directly  pertinent  to  the  subject  of  measurement  uncertainties  (a  subject 
considered  in  Appendix  II), 


2.  TIIK  TKMI'OH \I,->I’\TI\I.  HKSOI.I  TION  01  R\I»\R  MK  \SI  H  KM  KM'S  MVI1K  KITH 

pi  i, si:  or  muko  intiii.r  vhon 

The  signal  stability  that  is  achieved  by  video  integration  is  a  desirable  and 
essential  feature  of  weather-radar  measurements.  It  permits  the  average  signal 
levels  within  storms  to  be  determined  and  displayed  in  contoured  form  without  the 
presence  of  a  lot  of  fine-scale,  echo  detail,  which  otherwise  would  make  contouring 
difficult,  if  not  impossible. 
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It  must  be  emphasized,  however,  that  this  signal  stability  is  achieved  at  the 
cost  of  a  degradation  of  both  the  temporal  and  spatial  resolution  of  the  radar 
measurem ?nts  (relative  to  their  non-integrated  counterparts). 

The  fundamental  time  resolution  of  a  radar  with  hydrometeor  targets  is  the 
time  duration  of  the  backscattered  pulse,  which,  reference  Section  2.4.  1,  is 

At  -y..  (Al) 

where  h*  is  the  pulse  duration,  specified  in  time  units. 

The  fundamental  spatial  resolution  of  the  radar  is  that  of  the  pulse  volume, 
which,  from  Eq.  (28)  rewritten,  is 

V  =  >1  .  (A2) 

where  h  is  the  pulse  length  in  length  units.  The  relation  between  h  and  h'  is  given 
in  Eq.  (24). 

With  pulse  integration,  or  video  integration,  the  time  resolution  of  the 
measurements  becomes 


AtI  iS-* 


w  here  n^  is  the  number  of  radar  pulses  included  in  the  integration  and  prf  is  the 
pulse  repetition  frequency  of  the  radar. 

Ignoring  derivation  details,  the  spatial  resolution  of  measurements  made  with 
a  pulse  or  video  integrator  is 


'i  V  [l  f^rtan(^v]  • 


where  n^  is  the  number  of  range  cells  integrated  and  /i  is  the  antenna  scan  rate 
(for  example  ,i,  for  Rill  scan,  is  the  rate  of  change  of  the  elevation  angle, 
whereas,  ,i,  for  PPI  scan,  is  the  rate  of  change  of  the  azimuth  angle).  This 
equation  applies  only  to  a  parabolic  antenna. 

Resolution  values,  both  fundamental  and  with  integration,  are  shown  in 
Table  Al  and  were  computed  from  Eqs,  (A2)  through  (A4)  for  a  radar  having 
h'  0.5  jisec,  0  -  0.48°,  prf  060  pulses'sec  and  ,1  -  l°/sec.  (This  corresponds 
to  the  FPS-18  radar  of  the  SAMS  program  in  Hill  scan,  see  Section  4  and  Table  8.) 
It  was  also  assumed  that  the  integration  %vas  performed  for  256  pulses  and  over 
8  range  cells.  (These  are  typical,  normal  values  for  the  AFCRL  integrator.)  The 
resolution  values  of  the  table  are  shown  for  five  different  radar  ranges. 
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Table  Al.  Comparison  Example  of  the  Temporal  and  Spatial  Resolution  of  the 
Wallops  Island,  FPS-18  Radar  in  RHI  Scan,  With  and  Without  Video  Integration 
|  ref.  Table  3  and  Eqs,  (A2)  through  (A4)| 


Slant  Range 
r 

Fundamental  Radar 

Resolution 

Resolution  With  Video  Integration 

(for  n  250,  n  8  and  ,j  1°  sec) 

p  r 

<NM) 

(m) 

Time  Duration  of 
Back-Scattered 
Return  Pulse 

At 

(sec) 

Pulse 

Volume 

V 

(m  S 

Integration  Time 
AtI 
(sec) 

Integration  Volume 
VI 

(mS 

5 

9,  200 

0.25  ■  10*G 

0.35  X  10° 

0.207 

4.  78 

*-  io(i 

10 

18,  5.30 

0.25  <  10*G 

1.42  X  10° 

0.  207 

19.41 

•  ioG 

15 

27, 800 

0.25  <  10*G 

3.  19  \  10(> 

0.  207 

43.  59 

■  10° 

20 

37,  000 

0.25  x  10’6 

5.68  <  10*’ 

0.  207 

77.  02 

■  10G 

25 

40,  330 

0.25  >  10"° 

8.87  -  10f> 

0.  207 

121. 20 

■  10G 

It  is  seen  that  the  temporal  resolution  of  the  radar  measurements  is  degraded 
from  0.  25  ^sec  to  0.207  sec,  with  integration,  a  factor  of  approximately  10G,  and 
that  the  spatial  resolution  at  all  ranges  is  degraded  by  about  a  factor  of  1  i.  The 
specific  amounts  of  the  degradation  for  any  radar  and  integrator,  are  given  by  the 
ratios  At  At,  and  V  V,,  which  from  Eqs.  (Al)  through  (A4),  are 


At  JPprf 
At,  2np 


and 


V 


[1  +  ^tan<'^v] 


(A  5) 


(AO) 


independence  \m>  represent  \ti\ene>>oe  the  iimirdwetedr  sxmpi.es 

DBT  VINE-1)  II)  l’l  I.SE  OH  MOEO  INTEI.RVIIOV 


Assuming  for  the  moment  that  the  antenna  is  stationary  and  that  there  is  no 
wind,  for  any  single  range  cell  of  the  integrator,  the  "reshuffling  time"  of  the 
hydrometeors  required  to  give  non-autocorrelated,  independent  data  is  the  time  r 


*These  subjects  are  treated  in  brief,  summary  fashion  herein.  For  a 
detailed  discussion,  reference  Sirmans  and  Dovtak.1 

1.  Sirmans,  D„  and  Doviak,  R.J.  (1973)  Meteorological  Radar  Signal 
Intensity  Estimation.  NOAA  T.  M.  ERE  NSSL-04,  National  Severe  Storms  Tab. , 
Norman,  Oklahoma,  September  1973 
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specified  by  Eq.  (41)  which  is  based  on  the  work  of  Atlas.  “  Therefore,  if  the 
integration  period  is  Atj,  it  follows  that 


(A  7) 


where  n.  is  the  number  of  independent  samples.  Hence,  from  Eq.  (41)  and 
Eq.  (43), 


583  n 

ni  prf  \  ’ 


(A8) 


for  X  specified  in  centimeters.  It  also  follows  that  the  number  of  radar  pulses 
per  independent  sample  is 

n 

i  .  =  0.00171  prf  X  ,  <A9) 

from  Eq.  (A8). 

f  or  the  radar  example  cited  earlier,  which  reflected  the  SARIS  conditions 
of  operation  with  the  FPS-18  radar  and  AFCRL  integrator,  the  number  of  in¬ 
dependent  samples  with  a  stationary  antenna  and  no  wind  is  14.fi,  from  Eq.  (A8), 
for  np  25(1,  prf  OfiO  and  X  10.  7  cm,  and  the  number  of  pulses  per  independent 
sample  is  17.  fi  from  Eq.  (A9).  This  indicates  the  magnitude  of  the  n^  and  i  . 
values  for  a  single  range  cell. 

We  may  now  consider  the  situation  when  the  antenna  is  moving. 

If  the  radar  antenna  moves  with  a  scan  rate  >1,  new,  independent  samples  are 
acquired  at  any  given  range  distance  whenever  the  angular  travel  of  the  antenna, 
relative  to  some  starting  reference,  has  become  equal  to  the  beamwidth.  Hence, 
with  the  antenna  in  uniform  motion,  new,  independent  samples  are  acquired  at 
the  rate  ( '/fi  and,  during  the  integration  period  Atj,  the  number  of  new’  sampios 
acquired  is 

Mt. 

Ant  — ^  .  (A  10) 

rn 

The  acquisition  rate  is  i'./0  to  a  close  approximation,  ignoring  second  order 
terms  of  minor  contribution. 

2.  Atlas,  D.  (19(14)  Advances  in  radar  meteorology,  from  Advances  in 
Geophysics,  10:317-478,  Academic  Press,  New  York. 
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The  wind  movement  through  any  single  range  cell  will  also  introduce  new, 
independent  samples  of  hydrometeors.  The  rate  of  introduction  is  given 
approximately  by 


Rn 


4|vhI 

h  +  2r  0 


(All) 


where  )Vj|j  is  the  magnitude,  or  "speed",  of  the  horizontal  wind  vector.  This 
equation  assumes  that  the  wind  transports  the  hydrometeors  across  the  "average" 
dimension  of  the  radar  range  cell,  which  has  the  range-length  h/2  and  the  beam- 
cross -section-diameter  rf).  (Such  assumption  is  justified  because  the  number  of 
independent  samples  introduced  by  wind  translation  is  relatively  small  and  of 
second  order  contribution. ) 

During  the  integration  period  Atj,  the  number  of  new,  independent  samples 
acquired  by  wind  transport  is 


Atij  RTAtj  . 


(A12) 


The  total  number  of  independent  samples  per  single  range  cell,  in  the  case 
of  both  antenna  motion  and  wind  translation,  is  given  by 


ni  -  n.  +  A  n,  +  An. 
mw  m  Sv 


(A13) 


which,  from  Eqs,  (A3),  (AH)  and  (A10)  through  (A12),  becomes 
n,  At. 


mv 


T 


m+A,  liZiil 

X  0  h  +  2r  n 


(A14) 


The  number  of  radar  pulses  that  pass  through  each  of  the  independent  samples, 
for  the  moving  antenna,  wind  translation  case,  is 


m  n. 


(A  15) 


mw 


or,  from  Eqs.  (A3)  and  (A  14), 
_ prf 


m  585  +  JLfli 

X  0  h  +  2r0 


(Alt!) 


Ji! 
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Limits  are  imposed  on  the  number  of  pulses  per  independent  sample  per 
integration  period  that  can  be  obtained  by  pulse  or  video  integration.  If  we 
recognize  that  a  certain  minimum  number  of  pulses  per  sample  per  integration 
period  is  required  for  data  representativeness  this,  in  turn,  imposes  limits  on 
the  antenna  scan  rate  that  can  be  used  under  different  conditions  of  wind  transla¬ 
tion.  These  limitations  may  be  demonstrated  as  follows: 

If  it  is  stipulated  that 


(A17) 


where  i  ^  is  the  minimum  lumber  of  radar  pulses  per  sample  per  integration 
period  needed  for  representa.iveness,  then,  from  Eqs.  (A1G)  and  (A17),  the 
representativeness  requirement  for  any  given  radar  and  wind  condition  is. 


v  ■  < 

nun  = 


585  ^  fj  1  H 
~  IT"  H+Trfl' 


and  the  associated  limitation  on  the  antenna  scan  rate  is,  from  Eq.  (A  18),  solved 
for  (3, 


r  .  n  JEl£  585  _  4'vIjl 
1  min  X  h  +  2r0 


The  maximum  number  of  radar  pulses  per  independent  sample  per  integration 
period  that  can  be  obtained  with  a  given  radar  is  the  number  for  the  case  of  a 
stationary  antenna  with  no  wind.  This  maximum,  from  Eq.  (A1G),  is  specified  by 


.  Aprf 
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which  shows  that  a  change  in  the  maximum  can  only  be  effected  by  changing  either 
the  prf  or  the  wavelength  of  the  radar. It  is  also  seen,  by  reference  to  Eq,  (A1G), 
that  any  antenna  motion  or  wind  translation  of  the  hydrometeors  will  cause  i  ^  to 
have  values  smaller  than  this  maximum. 


■  One  method  of  increasing  this  maximum  value  is  by  "frequency  diversity", 
in  which  the  wavelength  is  changed  by  an  amount  exceeding  the  receiv'd'  bandwidth 
from  pulse  to  pulse  of  the  radar.  Reference  Section  4  for  additional  comments. 
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The  significance  of  the  above  equations  may  be  illustrated  by  citing  a  parti¬ 
cular  example.  Assume,  first,  that  at  least  15  radar  pulses  per  independent 
sample  are  required  for  representativeness,  such  that  r  .  =15,  Assume,  next, 

that  the  radar  of  concern  has  X  =  10.7  cm,  0  -  0.48°,  h  150  m  and  prf  9G0 
pulses  sec  (corresponding  to  the  FPS-18  radar  of  the  SAMS  program). 

For  this  radar,  Kq.  (A20)  shows  that  the  maximum  obtainable  number  of 
pulses  per  independent  sample  is  17.  G,  which  exceeds  the  number  ^  15 

that  we  assumed  were  required  for  representativeness. 

Suppose,  now,  that  the  antenna  is  moving  but  there  is  no  wind.  Then, 

Eq.  (A19),  for  i  .  15,  becomes 

nun 

,1  4.48°  sec  (A21) 

which  stipulates  that  the  antenna  cannot  be  scanned  more  rapidly  than  this,  if 
representative  data  are  to  be  acquired. 

Assume,  finally,  that  the  speed  of  the  horizontal  wind  at  the  altitude  of  a 
particular  range  cell,  located  r  -  5  km  from  the  radar  site,  is  30  m  sec. 

Equation  (A  19),  for  this  wind  situation,  specifies  that 

,1<  4.23°  sec,  (A 22) 

to  obtain  representative  data.  The  change  in  ,i  from  Eq.  (A21)  is  very  slight, 
emphasizing  the  fact  that  the  wind  translation  term  in  the  previous  equations  is  of 
second  order,  at  least  for  the  FPS-18  radar  of  our  example. 

The  discussion  thus  far  has  been  limited  to  the  case  of  a  single  range  cell. 

If  multiple  range  cells  are  included  in  the  integration,  then,  since  each  cell 
has  been  defined  and  set  equal  to  the  range  extent  of  the  pulse  volume,  the  data 
for  any  individual  cell  are  independent  of  the  data  for  adjacent  cells. 

Therefore,  the  total  number  of  independent  samples  included  in  the  pulse  or 
video  integration  is  the  multiple  of  the  number  of  cells,  n  ,  times  the  number  of 
independent  samples  per  sin  jle  cell,  or 

n.  ncnj  ,  (A23) 

lT  raw 


which,  from  Eqs.  JA3)  and  (A  14),  becomes 
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for  X  In  centimeters.  * 

The  other  equations,  Eqs.  (A15)  through  (A20),  regarding  the  number  of 
radar  pulses  per  independent  sample  per  integration  period  and  the  restrictions  on 
the  antenna  scan  rate,  remain  unchanged  for  the  case  of  integration  over  multiple 
range  cells. 

I.  \  \RIWI.K  Kyi  \  IION' 

The  importance  of  the  equations  just  derived  to  the  estimation  of  the  variance, 
and  uncertainties,  or  radar  measurements  made  with  pulse  or  video  integration 
will  now  be  discussed. 

Marshall  and  Hitchfield^  have  shown  that  the  distribution  of  the  averages  of 

n  independent  samples  of  a  parameter,  such  as  P  which  is  the  parameter 

measured  by  the  pulse  or  video  integrator,  reference  Section  3.2,  will  vary  with 

n  as  illustrated  in  Figure  Al.  The  dashed  and  dotted  lines  show  the  envelopes  of 

the  first  and  second  standard  deviation  of  the  P  averages  about  the  n  <*  value 

of  P~ ,  which  is  shown  bv  the  solid  horizontal  line.  It  is  seen  that  the  variation 
r’ _ 

of  the  p  ,  averages  decreases  with  n,  as  would  be  expected. 

If  the  standard  deviation  of  any  large,  unspecified  number  of  single 
independent  radar  samples  of  hydrometeors  is  ojT^,  ,  and  if  n  such  independent 
samples  are  group-averaged  by  a  pulse  or  video  integrator,  the  standard  devia¬ 
tion  of  the  integrated  values  of  F  will  vary  with  n  approximately  as 

w  op- 


Reference  Atlas"*  or  Sirmans  and  Doviak, 1  for  example,  where  y  is  a  constant 
of  proportionality  (having  values  close  to  unity)  that  depends  on  such  tilings  as  the 
type  of  radar  receiver  used,  the  particular  nature  of  the  integration  process, 


There  are  t.vo  points  we  wish  to  emphasize  by  writing  the  above  equation. 

The  first  is  that  the  standard  deviation,  or  other  statistical  measure  of  the 

_  _  -1/2 

variance  of  P  (or  I),  will  decrease  approximately  as  n  '  for  hydrometeor 


variance  of  P^  (or  I),  will  decrease  approximately  as  n  '  for  hydrometeor 
targets  in  atm  ospheric  clouds  or  stoims.  The  second  is  that  liq.  (A25)  provides 

Actually,  this  will  differ  somewhat,  frem  integrator  to  integrator, 
dependent  on  the  details  of  the  input  quantization  and  averaging  processes, 
reference  Sirmans  and  Doviak.  1 

3.  Marshall,  J.S.,  and  Hitchfield,  \V.  (1953)The  interpretation  of  the 
fluctuating  echo  for  randomly  distributed  scatterers,  Part  I,  Can.  J.  Phvs. 
31:902-994.  '  ‘ — 
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n  -  Number  of  Independent  Samples  per  Integrction  Group 


Figure  Al.  Standard  Deviation  of  the  Distribution  of  Averages  of  n  Independent 
Values  of  P[  t  The  horizontal  solid  line  shows  the  average  corresponding  to 

n  ■  ar  .  The  dashed  isolines  labeled  o—  show  the  envelope  or  the  first  standard 
_  r 

delation  of  the  P  ,  averages  for  grouped  independent  samples  containing  u 

samples  per  group,  for  n  values  ranging  from  1  to  100.  The  dotted  isolines 
show  the  envelope  of  the  second  standard  deviation.  The  two  open  circles  show 
the  particular  op-  values,  referenced  in  the  text,  which  pertain  to  the  P 


averages  for  single  independent  samples.  The  scale  values  of  the  ordinate  are 
arbitrary;  they  will  differ  from  one  averaging  or  integration  situation  to  another. 
Derived  from  Atlas  - 


a  means  of  estimating  the  uncertainties  (indeterminancies)  of  radar  integration 
data  that  were  acquired  for  known,  or  ascertainable,  values  of  n. 

To  be  more  specific  concerning  the  latter  point,  it  was  demonstrati  d 
previously  that  flip  number  of  independent  samples  for  pulse  or  video  integration 
with  a  moving  antenna  was  given  by  Kq.  (A24).  We  may  substitute  this  equation 


into  Kq.  (A2a),  replacing  n  In  n.  ,  to  obtain 

lT 
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This  s how s  that,  for  any  given  radar  having  a  particular  wavelength  X,  beam- 
width  0,  and  pulse  repetition  frequency  (,.rf),  the  standard  deviation  of  the  P 
(or  D  values  measured  with  a  pulse  or  video  integrator  will  decrease  (1)  with 
the  number  of  pulses  integrated,  that  is,  with  the  integration  period,  see  Eq.  (A3), 
(2)  with  the  number  of  range  cells  included  in  the  averaging  and  (3)  with  the  scan 
rate  of  the  antenna  and  the  wind -transport  function  (within  the  scan  rate  limita¬ 
tions  of  Eq.  (A19)| . 

If  a  given  radar  has  a  selectable  prf,  it  is  likewise  apparent  from  Eqs  (A3) 
and  (A2f>),  that  since  the  integration  period  Atj  n  /  prf,  then  the  largest  select¬ 
able  prf— when  combined  with  a  long  integration  period  (which  is  a  matter  of 
integrator  design  and  or  operational  choice)— will  tend  to  minimize  the  ap^ 
deviations,  assuming  that  the  values  of  the  other  parameters  in  Eq.  (A2G)  are 
held  constant.  Conversely,  the  selection  of  a  small  prf  combined  with  a  short 
Atj  will  tend  to  maximize  the  opr  values.  These  statements,  although  true  for 
the  restricted,  normal  range  of  integrator  operation,  do  have  limitations  and 
should  not  be  assum’d  to  hold  generally. 

Equation  (A28)  can  be  employed  to  estimate  the  uncertainties,  or  indeter¬ 
minancies,  of  integrator  measurements  acquired  from  hydrometeor  targets. 

For  example,  if  opr  can  be  evaluated  for  single  independent  samples,  and  if 

l,  is  known,  reasonable  estimates  of  op  can  be  acquired  from  Eq.  (A20). 

Alternately,  if  op  can  itself  be  established  for  a  particular  integration  situation 

1  r 

of  known  n,  then  Eq.  (A2<i)  can  also  be  solved.  These  possibilities  are  merely 
suggested  here  and  will  not  be  elaborated  (see  Appendix  11  for  additional  dis¬ 
cussion).  fine  of  the  problems  involved  is  that  the  atmospheric  volumes  used 
for  such  estimates  must  contain  hydrometeors  that  are  relatively  homogeneous 

in  type,  quantity  and  size  distribution.  Otherwise,  the  cm  or  ott-  values  w  ill 

rl  1  r 

reflect  real  atmospheric  variability,  rather  than  just  t lie  uncertainty  values 
desired. 

The  work  of  this  appendix  is  summarized  in  Section  2.  7.4  of  the  main  text. 
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Appendix  B 

Parameters  Contributing  to  the  Uncertainties  of  Radar 
Measurement  and  Proposed  Uncertainty  Analysis 


The  parameters  that  contribute  to  the  uncertainties  of  radar  measurement 
may  be  identified  by  writing  the  equations  of  relationship  between  the  basic- 
measured  parameter  of  the  radar,  P  and  the  desired  end-product  parameters, 
reflectivity  factors  Z^  and  Zj,  for  water  and  ice  hydrometeors. 

As  noted  in  Sections  2.7.  1  and  2.7.2,  the  AFt'HF  integrators  employed  with 
the  KARF  radars  basically  measure  the  average  power  of  the  baekscattcr  signals 
received  from  the  hydrometeor  particles  that  exist  within  the  integration  volume 
of  the  radar  beam  during  the  particular  period  of  integration.  This  average 
icccived  power,  P^,  is  related  to  the  true  average  received  power,  1*,.^,  as 
specified  by  Fq,  (4‘D.  The  integration  signal,  I  ,  produced  by  the  video 
integrator,  is  related  to  P  as 

1  e'Pj.r2  .  (Ill) 

See  Fq.  51.  The  volume  reflectivity,  »j,  is  dependent  on  I  in  the  manner 


See  Fq.  (40).  Finally,  the  end-product  parameters  of  the  SAMS  radar  measure¬ 
ments  program,  the  reflectivity  factors  Z^  for  water  hvdrometeors  and  Zj  for 
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ice  hydrometeors,  are  related  to  n  as  specified  respectively  by  laps,  (5  1)  and 
(f.4). 

By  the  successive  substitution  of  Kq.  (43)  into  Kq.  (HI),  Kq.  (HI)  into  1  q.  (B2) 
and  Kq.  (B2)  into  Kq.  (57),  first,  and  Kq.  (B2)  into  Kq.  (04)  second,  it  is  readnv 
demonstrated  tliat 


for  water  liydrometeors,  and 


(B4) 


l  •> 

for  ice  liydrometeors.  The  multiplying  factor  10  "  has  been  incorporated  in  these 
equations  so  that  Z^.  and  Zj  will  be  given  in  the  conventional  units  of  mm'1  ni  , 
for  A  in  cm  and  r  in  N'M. 

Kquations  (B3)  and  (B4)  show  the  relationships  between  Z^.  and  Prt,  and  Zj 
and  P  ,  which  contain  all  of  the  factors  of  dependence.  These  are  the  factors 
that  are  individually  subject  to  uncertainty  and  that,  in  their  total  uncertainty 
effects,  will  cause  uncertainty  of  the  Z^.  and  Zj  values  for  the  SAMS  missile 
trajectories.  Kach  of  these  factors  will  be  discussed,  in  turn,  in  the  following 
paragraphs. 

The  factors  r,  A  and  c  are  subject  to  only  minor  uncertainties  and  can  be 
regarded  as  essentially  devoid  of  error  contribution.  The  radar  range  r,  is 
determinable  to  at  least  a  tenth  of  a  microsecond  in  time  or  15  m  in  range  dis¬ 
tance.  The  radar  wavelength  A,  during  radar  operation,  does  not  vary  by  more 
than  about  1  or  2  percent.  When  the  video  integrator  is  coupled  »o  the  radar 
receiver,  the  coupling  constant,  c,  can  be  set  to  within  at  least  a  tenth  of  a 
decibel  accuracy,  which  is  about  t  1  percent.  These  uncertainties  are  definitely 
of  second  order  compared  to  the  other  uncertainties  that  enter  Kqs.  (B3)  and  (B4). 
The  numerical  values  of  the  averaging  factor,  ■>,  used  in  the  computational 

equations  of  the  main  text,  stem  from  the  work  of  Kodaira,  *  Austin  and 
'> 

Schaffner,  “  and  perhaps  others.  The  values  depend  on  assumptions  regarding 


1.  Kodaira,  X.  (1000)  The  Characteristics  of  the  Averaged  Kcho  Intensity 
Received  by  the  Logarithmic  1. 1'.  Amplifier.  Proceedings  8th  Weather  Radar 
Conference,  Am.  Meteorol  Soc. ,  Boston,  Massachusetts,  pp  255-201. 


Austin,  P.  M. ,  and  Schaffner,  M.R.  (1070)  Computations  and  Kxperi- 
telcvant  to  Digital  Processing  of  Weather  Radar  Kchocs.  Proceedings 


14th  Weather  Radar  Conference,  17-20  November  1070,  Tucson,  Arizona. 
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tl:e  probability  that  the  average  of  n  independent  samples  of  hydrometeors  detected 

_ o 

by  radar,  which  liave  a  mean-square -voltage  return  F“  (or  average  received 

power T*  ,)  will  exceed  a  certain  threshold  level,  which  lies  a  dll  amount,  y  , 

r  _  •)  _  1 

below  the  true  average  value  Et“  (or  Prj).  Thus,  y,  or  Lyj,  depends  on  the 

number  of  independent  samples  integrated  and  also  depends  on  the  particular 
types,  and  instrumental  characteristics,  of  the  radar  receivers  and  video 
integrators  employed.  The  estimation  of  the  assumptive  errors  involved  in  the 
c  hoice  of  a  specific  value  of  y  is,  therefore,  seen  to  involve  considerable 
complexity.  The  errors  are  probably  of  the  order  of  £0.2  to  tl.O  dB  and  it  is 
apparent  that  these  will  be  "bias  errors",  that  will  act  to  shift  all  of  the  Z^.  or 
Zj  trajectory  values  either  upward  or  downward  by  some  la  to  ±23  percent  (about 
an  assumed  typical  value  of  y  0.  537). 

The  hydrometeor  shape  factor  k  (see  Section  2.8)  was  assumed  equal  to  unity 
'n  the  present  report,  both  for  water  hydrometeors  and  ice  hydrometeors.  Tor 
water,  that  is  water  cloud  or  rain,  this  means  that  the  water  drops  are  assumed 
to  have  spherical  shape.  Studies  liave  shown  that  cloud  and  raindrops  up  to  about 
1-mm  diameter  size  are  spherical  within  close  tolerance.  From  the  1  mm  size 
up  to  the  "drop  breakup  size'  of  approximately  5  mm,  however,  the  falling  drops 
tend  to  depart  more  and  more  from  sphe  deity  and  to  develop  oscillatory 
characteristics  in  which  their  dimensions,  vertical  versus  horizontal,  become 
variable  with  fall  distance.  On  the  average,  though,  even  for  these  larger  droplets, 
the  k  1. 0  assumption  would  appear  to  hold  quite  well.  For  example,  Newell  and 
Cieotis,  *  from  backscatter  measurements  made  with  horizontally  and  vertically 
polarized  radiation,  found  no  disccrnable  differences  attributable  to  shape  factor 
in  any  hydrometeor  region  other  than  the  melting  zone.  (There  was  a  1  db  differ¬ 
ence  in  this  zone  indicating  preferential  horizontal  orientation. )  Other  information 
concerning  this  subject  undoubtedly  exists  in  the  literature. 

Likewise,  the  values  of  the  refractive-index -factor,  ]  K  for  water  hydro - 

meteors,  are  known  with  considerable  precision.  The  values  vary  with  the  radar 

q 

wavelength,  as  discussed  by  Gunn  and  Fast,  but  there  is  little  question  about 
their  accuracy  and  theoretical  justification. 

For  all  intents  and  purposes,  then,  the  k  and  !  K !  ^  parameters  for  water 
hydrometeors  may  be  regarded  as  devoid  of  important  error  contribution  to  Z^.. 


3.  Newell,  R.F. ,  and  Goetis,  S.  G.  (1055)  Meteorological  Measurements  With 
a  Radar  Provided  With  Variable  Polarization.  M.I.  T.  Dept,  of  Meteorol.  Weather 
(tadar  lies. ,  Tech,  Note  No.  8, 

4,  Gunn,  K.  L.S.,  and  Fast,  T.D.ll.  (1954)  The  microwave  properties  of 

pre<  ipitation  particles,  Quart.  .1.  Bov.  Met.  Soc.  60:522.(427-8). 
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In  the  rase  of  ice  hydrometeors,  particularly  ice  crystals,  it  is  not  a-priori 
apparent  that  the  shape  factor,  k,  and  the  density-modified-refractive-index- 

.  |  9 

factor,  K' ,  ",  are  devoid  of  uncertainties.  The  particle  shapes,  except  for  some 
grauple  particles  and  hail,  are  mostly  non-spherical.  The  ice  crystals  have 
various  needleform,  columnar,  platelet,  dendritic  and  irregular  shapes,  and  the 
snow  particles  are  composed  of  aggregates  of  such  crystals.  The  k  1.0 
assumption,  as  noted  in  Section  2.  0,  is  contended  to  be  justified  for  these  non- 
spherical  particles  on  the  basis  that  they,  while  falling  within  the  radar  pulse 
volume,  or  integration  volume,  are  probabilistically  likely  to  have  random 
orientations  in  space.  This  assumption,  however,  at  least  for  the  case  of  a  stable 
atmosphere  with  little  atmospheric  turbulence  such  as  would  lie  expected  to  prevail 
in  the  Wallops  storms  desired  for  SAMS  missile  launch  purposes,  is  in  conflict 
with  various  observational  and  laboratory  findings  that  ice  crystals  and  single¬ 
crystal-snow  particles,  except  for  the  irregular  and  dendritic  forms,  have 
preferential  fall  orientations.  Thus,  it  would  appear  that,  although  the  k  1.0 
assumption  should  be  reasonably  valid  for  all  types  of  ice  particles  under  turbulent 
conditions,  it  might  be  questionable  for  particles  with  preferred  fall-orientations 
under  stable,  relatively  nonturbulent  conditions.  The  errors  would  conceivably  be 
ua  large  as  a  factor  of  two,  in  the  author's  opinion.  •  Detailed  data  concerning  the 
horizontally -polarized  versus  vertically -|jolarized  signal  returns  from  such  ice 
crystals  under  stable  conditions  should  be  acquired  to  provide  the  necessary, 
error -assessment  information. 

The  density-modified-refractive-index-factor  ^  K'  |  “  for  ice  hydrometeors, 
which  was  assigned  ‘.he  value  0.20!)  herein  and  was  based  on  the  work  of  Mason'’ 
(see  Section  2,8),  would  seem  subject  to  relatively  minor  uncertainty  perhaps  of 
the  order  of  t4  percent.  The  cited  value  was  obtained  under  the  assumption  that 
the  ratio  of  the  refractive  index  of  ice,  Kj,  to  the  density  of  ice,  p^  would  be 
constant  for  ice  of  all  densities  (that  is,  for  ice  containing  different  amounts  of 
entrained  air  bubbles).  The  assumption,  for  single-|>article,  ice  crystals,  is 

7 

quite  valid  (see  Atlas  et  al  )  since  the  values  of  Kj,  pj  only  vary  from  0.  4G  to  0.  r>0 

Rased  on  the  work  of  Atlas  et  al,  l>  concerning  scattering  from  non-spherical 
particles  in  the  Rayleigh  region,  also  based  on  the  comments  of  Mason.  ” 

f>.  Atlas,  I).,  Kerker,  M. ,  and  Hitschfcld,  W  .  (1933)  Scattering  and 
attenuation  bv  non-spherical  atmospheric  particles,  J.  Atmos.  Terr.  Phvs.  3: 

108.  (400-13).  ~ 

0.  Mason,  15.  J.  (1071)  The  Physics  of  Clouds,  Second  Edition,  Clarendon 
Press,  Oxford,  England. 

7.  Atlas,  D. ,  Cunningham,  R.M.,  Donaldson,  R.J.,  Jr.,  Kantor,  G. ,  and 
Newman,  P.  (190.'))  Some  aspects  of  electromagnetic  wave  propagation.  Handbook 
of  Geophysics  and  Soacc  Environments.  9-1  through  9-20,  AFCRI.,  Bedford, 
Massachusetts. 
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for  p.  values  spanning  the  normally -observed,  ice-densitv  range  from  0.22  to 

1  -3 

0.  97  gm  cm  .  The  approximation  is  also  valid  for  Hayleigh  scattering  from 
snow,  since  snow  may  be  considered  to  be  a  mixture  of  the  two  dielectrics,  ice 
and  air.  The  composite  K  p  ratio  for  such  mixture  is  given  by  the  mass -weighted 
sum  of  each  substance.  Hence,  since  K  p  for  air  is  negligible  compared  to  that 
for  ice,  the  assumptive  uncertainties  for  snow  particles  are  not  appreciably 
larger,  theoretically,  than  for  single-particle,  ice  crystals. 

The  uncertainty  n  I  system  calibration  constant,  k,  were  discussed 
previously  in  Section  3.  5.  The  difficulties  of  assessing  the  uncertainties  were 
pointed  out  and  methods  for  acquiring  appropriate  uncertainty  information  were 
suggested.  As  a  first,  gross  estimate,  it  was  contended  that  the  uncertainties 
were  probably  of  the  order  of  tO.  dll  to  t3.  5  dll  (■*■  t  12  percent  to  t  120  percent).  * 
This  contention  is  in  reasonable  accord  with  the  Stout  and  Mueller  estimation 
that  under  ideal,  i  arcful  conditions  of  radar  calibration,  the  accuracy  is  about 
t3.0  dB.  Mason*’  cites  accuracies  of  approximately  ±2  dll.  It  should  be  noted 
that  these  calibration  uncertain  lies  are  of  "bias  type",  which,  for  any  given  day 
of  radar  operation,  will  cause  a'l  of  the  or  Zj  values  for  the  missile 
trajectories  to  be  either  too  large  or  too  small. 

To  this  point,  we  have  considered  the  uncertainty  status  of  all  factors 
entering  Eqs.  (B3)  and  (B4)  except  that  of  the  basic  parameter  of  measurement, 
Prt.  It  will  be  convenient  at  this  time  to  rewrite  Eqs.  (B3)  and  (B4)  and  make  a 
clear  distinction  between  those  parameters  which  are  not  subject  to  important 
uncertainty  and  those  that  are.  Thus,  10 q.  (B3)  may  be  written  as 


4  2  l‘> 
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where  7,  k  and  P  are  the  uncertain  parameters,  and  Eq,  (B4)  may  be  written 

rt 

as 
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where  k  is  an  additional  factor  of  uncertainty  (for  the  particular  case  of  ice  crystals 
crystals  having  preferential  fall  orientations). 


Reference  footnote  on  page  102,  concerning  percentage  uncertainties. 

8.  Stout,  G,  E. ,  and  Mueller,  F.  A.  (1968)  Survey  of  relationships  between 
rainfall  rate  and  radai  reflectivity  in  the  measurement  of  precipitation,  J.  Appl. 
Meteorol.  7:465-74. 
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The  total  variation  in  Z^.  resulting  from  component  variation  of  ■>,  k  and 
is 
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while  that  for  Zj,  for  the  added  component  variation  of  k,  is 
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When  the  partial  derivatives  of  the  above  equations  are  evaluated,  from  Eqs. 
(1)5)  and  (Bfi),  we  obtain 
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The  derivatives  dy,  ck,  dk  and  dPrt  of  these  equations  may  be  dentified  and 
associated  with  the  uncertainty  variations  of  the  parameters  y,  k,  k  and  Prt  . 

Their  total  contributions  to  the  uncertainties  of  the  reflectivity  factors  for  water 
and  ice  hydrometeors,  that  is  dZ^.  anti  dZj,  are  then  specified  by  Eqs.  (119) 
and  (1)10), 

\\  e  turn,  now,  to  a  discussion  of  the  uncertainties,  or  indeterminancies,  of 

the  true  received  power,  P  , 

rt 

An  equation  was  presented  in  Appendix  A,  Eq.  (A25),  which  specified  how 
the  standard  deviation  of  the  integrated  received  power  would  vary  (approximately) 
with  the  number  of  independent  hydrometeor  samples  included  in  the  integration. 
Another  equation,  Eq.  A2fi,  was  also  written  which  ue^cribed  how  the  standard 
deviation  of  the  received  power  would  vary  with  the  particular  characteristics  of 
a  given  radar  and  with  the  integration  and  wind-translation  parameters. 

l  or  the  purposes  of  the  uncertainty  analyses  of  this  appendix,  we  may  assume 
that  the  uncertainty  parameter  dP  ,  of  Eqs.  (119)  and  (1110),  is  given  approximately 

by  1 


(Bll) 
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where  ap  is  as  specified  bj  Fq.  (A2C).  This  assumption  permits  Eqs.  (B9) 


and  (BIO)  to  be  written  as 
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where  a[7^  ,  as  defined  in  Appendix  A,  is  the  standard  deviation  of  the  received 
power  for  Single  independent  samples. 

These  are  the  final  uncertainty  equations  for  dZ^.  and  dZj.  * 

To  illustrate  the  nature  of  the  uncertainty  analyses  permitted  b>  these 
equations,  the  following  example  is  presented. 

Assume  that  we  wish  to  assess  the  uncertainties  of  the  Z,^  and  Zj  values  for 
the  missile  trajectories  that  were  determined  by  the  FPS-18  radar  and  integrator 
of  the  1072-7.1  season  of  SAMS  operations.  For  this  radar  and  season, 

X  10.7  cm,  fl  -  0. 48°,  h  150  m,  ->  0.  537  and  c  8.42  <  1011mw'1cm  1 

(N’M)*2,  see  Tables  3  and  4.  The  radar  calibration  constant,  «  ,  for  the  season 

17 

had  an  average  value  of  about  8.0-10  ,  see  Table  0.  The  integration  para¬ 

meters  assumed  were  n^  250,  prf  000,  |i  1°  sec,  and  u  1.0.  The  number 
of  range  cells  integrated,  nc,  is  left  unspecified  at  the  moment,  for  uncertainty 
investigation. 

It  should  be  noted  that,  if  frequency  diversity  is  used  in  the  video  integra¬ 
tion,  as  planned  for  the  Spandar  radar  in  the  1973-74  season,  the  number  of 
independent  samples  integrated  will  be  increased  considerably.  Such  change  will 
affect  the  last  terms  of  Fqs.  (1112)  and  (B13)  and  the  modifications  indicated  in 
Section  3.7.4  will  have  to  be  accomplished  to  obtain  new,  appropriate,  uncertainty 
equations. 


If,  in  addition  to  the  values  of  these  cited  parameters,  we  presume  that  the 

'  — ^  I  __  .  | 

normal,  average  wind  speed  in  the  W  allops  storms  is  I  VH>  15  m  sec  , 
the  uncertainty  equation  for  water  hydrometeors,  Eq.  (B 12),  for  which  k  1.0 
and  *k!2  0.93,  reduces  to 
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\\  rt  t 


0.537  Op 


0.2G7  n  f  5li.7  +  rr^ 

cl  loO  +  11). 


for  r  specified  in  km  and  in  mw. 

For  snow  hydrometeors)  or  for  ice  crystals  that  do  not  have  preferred  fall 
orientation,  for  which  k  1.0,  |  K'  I  “  0.209,  and  dk  ~  0,  the  uncertainty 

Eq.  (1113)  becomes 
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For  that  special  class  of  ice  crystals  which  do  have  preferred  fall  orienta¬ 
tion,  that  is,  for  which  dk  *  0,  Eq.  (1113)  reduces  to 

2 
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The  individual  term  contributions  to  the  total  uncertr inties,  d Z^  and  dZj, 
are  assessed  in  Tables  HI  through  114  for  four  hydrometeor  types:  (1)  rain, 

(2)  snow,  (3)  ice  crystals  without  pieferential  fall  orientation,  and  (4)  ice 
crystals  with  preferential  fall  orientation. 


Table  B3.  Estimated  Term  Contributions  and  Total  Uncertainties  of  the  Radar  Reflectivity  Factor  for  Ice  Crystals 


The  equations  are  evaluated  for  two  radar  ranges,  r  5  NM  and  r  10  NM, 
which  correspond  approximately  to  the  mid-range  and  extreme  range  of  the  missile- 
trajectory  measurements  of  the  1970-73  seasons.  Two  values  of  the  integration 
signal  level,  ^  are  specified  in  each  table  which  span  the  normally-observed  .  1 . 
values  of  the  particular  hydrometeor  type  in  the  Wallops  storms.  The  correspond¬ 
ing  Prt  and  Z^„  or  Zj,  values  are  indicated  which  were  computed  from  Eq.  (51) 
and  Eqs.  (B3)  and  (B4). 

The  assumptions  about  the  mean  values  and  uncertainties  of  the  independent 
parameters  of  Eqs.  (B14)  through  (B16)  are  consistent  with  the  previous  text  dis¬ 
cussion  of  this  appendix.  Specifically,  it  was  assumed: 

(1)  That  the  averaging  factor,  y,  is  uncertain  to  ±0.5  dB,  giving  dy  ±0.  062 
about  an  average  y  of  0.  537. 


d« 


(2)  That  the  calibration  constant,  k,  is  uncertain  to  ±2.0  dB,  giving 
“  3.  8  '<•  1017  about  an  average  k  of  8.  0  <  lO17. 


(3)  That  the  particle  shape  factor,  k,  is  devoid  of  uncertainty  for  rain,  snow 
and  ice  crystals  without  preferential  fall  orientation;  but  that  it  is  uncertain  to 
±2.0  dB  for  ice  crystals  with  preferential  fall  orientation  under  stable  atmospheric 
conditions,  giving  dk  “S'  ±0.48  about  an  average  k  of  unity. 

(4)  That  the  standard  deviation  of  the  integrated  received  power  for  single 
independent  samples,  opr  ,  is  uncertain  to  ±10  dB,  as  the  first  assumption 

indicated  in  the  tables,  and  is  uncertain  to  ±5  dB,  as  the  second  assumption 

indicated  to  the  tables;  the  uncertainties  being  with  respect  to  the  P  values 

1 1 

shown  in  the  third  columns. 

The  table  computations  were  also  performed  for  two  different  values  of  range¬ 
cell  integration,  for  single  range  cells,  n^  ■  1,  and  for  eight  cell  integration, 
n^  =  8.  Single  range-cell  integration  was  used  with  the  FPS-18  radar  in  the 
1970-71  season;  eight  range-cell  integration  was  used  in  the  1971-73  seasons. 

The  estimated  total  uncertainties,  dZ^v  (or  dZj)  of  the  next  to  the  last  columns 
of  Tables  B1  through  B4  w  ere  computed  as  the  probable  error  (P.  E.)  summation  of 
the  separate  term  contributions.  Actually,  for  any  given  period  and  day  of  SAMS 
radar  measurements  at  Wallops  Island,  the  uncertainty  factors  d>  and  d«  will 
impose  bias  errors  on  the  measurements.  But,  since  we  cannot  assume  a-priori 
knowledge  of  the  signs  of  these  factors  in  the  table-estimations  herein,  they  are 
treated  as  random  uncertainties  and  summed  in  the  P.E.  fashion. 

The  table  results  reveal  that  the  percent  uncertainties  of  Z^  (or  Zj)  vary 
from  about  ±51  percent  to  ±136  percent,  for  rain,  snow  and  ice  crystals  without 
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preferential  fall  orientation,  and  from  about  ±70  percent  to  ±144  percent,  for  ice 
crystals  with  presumed,  preferential  fall  orientation. 

The  percent  uncertainties  for  nc  *  8  are  seen  to  be  some  60  percent  to  70 
percent  smaller  ti  in  those  for  nc  =  1,  in  the  case  of  op  =  ±10dB,  and  some 
10  to  11  percent  smaller,  in  the  case  of  o pr  «  ±5dB. 

It  should  be  reiterated  that  Tables  Bi  through  B4  are  presented  primarily  to 
illustrate  the  techniques  of  uncertainty  analyses  based  on  Fqs.  (B12)  and  (B13). 
The  particular  table  values  are  merely  first,  gross  estimates  of  possibility.  For 
any  refined  uncertainty-assessment,  specific  data  will  have  to  be  acquired  for 
each  of  the  Wallops  storms,  concerning  the  values  of  the  independent  parameters 
that  enter  Fqs.  (B12)  and  (B13). 


Note  \ddrd  h olio** in^  Munu>i-ri|it  t  .iHnpli-linn 

The  careful  reader  will  note  that  the  P  values  of  Tables  Ill  through  B4  differ 
somewhat  from  those  specified  by  the  equation 


F  —LT.  (B17) 

rt  cy  r 

which  stems  from  the  combination  of  Eq.  (43)  and  Kq.  (Bl).  The  P  values  of  the 
tables  should  have  been  computed  from  the  table  X  values  using  this  above  equa¬ 
tion  (w  ith  appropriate  dB  to  absolute  conversion).  In  actual  fact,  however,  they 
were  inadvertently  computed  from  the  relationship 

f*  =  .  (13 1 8) 

rt  o  r 

which  neglects  the  y  factor  of  Eq.  (B17), 

The  effect  of  this  error  on  the  final  uncertainty  values  of  the  tables  is  rather 
slight  and  inconsequential,  particularly  in  view  of  the  Tact  that  the  table  uncer¬ 
tainties  are  merely  gross  estimates. 


Although  it  is  convenient  to  cite  percentage  uncertainties  for  qualitative 
comparison  purposes,  it  should  be  noted  that  large  uncertainties,  such  as  the  above, 
which  are  expressed  in  percent,  are  highly  non -symmetric  in  their  positive  versus 
their  negative  values.  For  example,  the  positive  numerical  value  ran  exceed  100 
but  the  negative  num?rical  value  cannot.  The  meaningful  uncertainty  values  of 
Tables  Bl  through  B4  are  the  d/.  values  of  the  next  to  the  last  columns,  the 
percent  uncertainty  values  of  the  last  columns  are  merely  indicative  and  must  lie 
interpreted  with  care. 


102 


References 


Kodaira,  N.  (1960)  The  Characteristics  of  the  Averaged  Echo  Intensity 
Received  by  the  Logarithmic  l,F.  Amplifier.  Proceedings  8tl  Weather 
Radar  Conference,  Am.  Meteorol.  Soc.,  Boston,  Massachusetts, 
pp  255-261. 

Austin,  P.M.,  and  Schaffner,  M.  R.  (1970)  Computations  and  Experiments 
Relevant  to  Digital  Processing  of  Weather  Radar  Echoes.  Proceedings 
14th  Weather  Radar  Conference,  17-20  November  1970,  Tucson,  Arizona. 

Newell,  R.E, ,  and  Geotis,  S.  G.  (1955)  Meteorological  Measurements  With 
a  Radar  Provided  W  ith  Variable  Polarization,  M.  I.  T.  Dept,  of  Meteorol. 
Weather  Radar  Res. ,  Tech.  Note  No.  8. 

Gunn,  K.  L.  S. ,  and  blast,  T.D.R.  (1954)  The  microwave  properties  of 
precipitation  particles.  Quart.  J.  Roy,  Met,  Soc.  80:522.  (427-8). 

Atlas,  I).,  Kcrker,  M. ,  and  Hitschfeld,  W.  ( 1958)  Scattering  and  attenuation 
by  non-spherical  atmospheric  particles,  J.  Atmos.  Terr.  Phys.  3:108. 
(4C6-13). 

Mason,  B.J.  (1971)  The  Physics  of  Clouds.  Second  Edition,  Clarendon 
Press,  Oxford,  England. 

Atlas,  D.,  Cunningham,  R.M.,  Donaldson,  R.J.,  Jr.,  Kantor,  G. ,  and 
Newman,  P.  (1965)  Some  aspects  of  electromagnetic  wave  propagation, 
Handbook  of  Geophysics  and  Space  Environments.  9-1  through  9-26, 
AFCRL,  Bedford,  Massachusetts. 

Stout,  G.  E. ,  and  Mueller,  F.  A.  (1968)  Survey  of  relationships  between 
rainfall  rate  and  radar  reflectivity  in  the  measurement  of  precipitation, 

J.  Appl,  Meteorol.  7:465-74. 


103 


List  of  Symbols 


\rubir  Symbols 


a 


A 


c 


C 


C 


r 


D 

De 

DI 

dB 

dBm 


decima1  attenuation,  per  unit  path  length 
attenuation  factor 
effective  antenna  area 

actual  cross-sectional  area  of  parabolic  antenna 
sampling  area  of  disdrometer  instrument 
geometric  cross-section  of  a  scattering  target 
c  oupling  constant 

calibrated  constant  of  the  radar  equation  pertaining  to  hydrometeor 
targets  with  video  Integration 

constant  of  reference  performance  for  the  radar  equation  pertaining 
to  hydrometeor  targets  with  video  integration 

diameter  of  water  droplets 

equivalent,  melted  diameter  of  ice  crystal  or  snow  particles 
physical  measure  of  the  actual  size  of  ice  crystals  or  snow  particles 
decibel 

decibels  below  t  milliwatt,  of  reference  power 
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—2 

E  mean  square  voltage  amplitude 

Gf  gain  of  a  receiving  antenna 

Gt  gain  of  a  transmitting  antenna 

h  pulse  length,  in  length  units 

h'  pulse  duration,  or  "pulse  length"  in  time  units 

I  integration  signal,  from  video  integrator 

1 ,  10  log  I,  thai  is,  decibel  value  of  1 

I  *  reference  integration  signal,  of  threshold  setting 

k  shape  factor,  for  hydrometeor  particles 

k'  microwave  absorption  coefficient  for  water 

K  antenna  calibration  constant,  of  "pattern  function"  equation 

Kq  antenna  calibration  constant,  of  "effective  area"  equation 

1  k]2  refractive  index  factor,  for  water  hydrometeors 

|  K'  |  2  density -normalized  refractive  index  factor,  for  ice  hydrometeors 

kr  sphere  calibration  constant  of  radar  reference  performance 

kg  sphere  calibration  constant  of  actual  radar  performance 

kg  10  log  kg,  that  is,  decibel  value  of  kg 

L  radar  system-loss-factor,  as  specifically  defined  in  the  text 

IV1  precipitation  liquid  water  content 

n  number  of  independent  samples,  general 

nc  number  of  range  cells  included  in  pulse  or  video  integration 

n.  number  of  independent  samples,  in  the  case  of  a  stationary  antenna 

and  no  wind 

n.  number  of  independent  samples,  with  antenna  motion  and  wind 

mw 

n.  total  number  of  independent  samples  included  in  the  pulse  or  video 

XT  integration 

np  number  of  radar  pulses  per  integration  period 

N  number  of  hydrometeor  particles  per  unit  volume 

Nj.  number  of  hydrometeor  particles  per  unit  area  per  unit  time 

NM  nautical  miles 
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Pr. 


nun 


‘  t 
prf 

R 

Rq, 

r 

r 
S 
S 

t 

<lt 

At 

V 

U 


H 


I 


V 

V, 


II1 


7 

7, 

zvv 


D 


'I 


received  power 

pulse-integrated  received-power 
true  pulse-integrated  received-power 

minimum,  integrated -received -power,  of  threshold  setting  for  the 
video  integrator 

transmitted  power 

pulse  repetition  frequency 

rainfall  rate 

wind  transport  rate,  of  hydrometeors  into  radar  range  cell 

radar  range 

correlation  coefficient 

power  density  in  space,  or  field  strength 

beam  shape  factor,  of  Probert -Jones 

time 

fundamental  time  resolution  of  the  radar 
time  resolution  with  pulse  or  video  integration 

units  conversion  factor  for  radar  equation  for  hydrometeor  targets 
units  conversion  factor  for  radar  equation  for  single  targets 
velocity  of  microwave  propagation  (3.0  7  1010  cm  sec"1) 
pulse  volume 

integration  volume,  for  pulse  or  video  integration 
magnitude,  or  speed,  of  the  horizontal  wind  vector 
circuit  impedance 
radar  reflectivity  factor,  general 

radar  reflectivity  factor  for  water  hydrometeors,  determined  from 
radar  measurements 

radar  reflectivity  factor  for  water  hydrometeors,  determined  from 
surface  disdrometer  measurements 

radar  reflectivity  factor  for  ice  hydrometeors 
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(»re**k  Smbols 


y 


e 

n 

o 

o 

K 


V 


t' 


o 


r 


0 

Id 

n 


wavelength  dependent  parameter  entering  the  calibration  equation 
of  the  disdrometer  method 

angular  scan  rate  of  the  antenna 

averaging  factor,  relating  mean  square  voltage  amplitude,  measured 
by  pulse  or  video  integrator,  to  average  received  power 

dielectric  constant,  for  hydrometeor  particles 

volume  reflectivity 

average  volume  reflectivity,  as  measured  with  a  video  integrator 
antenna  beamwidth,  for  parabolic  antenna 

radar  calibration  constant,  for  disdrometer  method  of  calibration 
10  log  k  ,  that  is,  decibel  value  of  k 

the  mean  value  of  determined  from  disdrometer  calibration 
radar  wavelength 

number  of  radar  pulses  per  independent  sample  per  integration 
period 

density  of  ice 

reflectivity  of  scattering  particles  for  microwaves 
density  of  water 

back  scattering  cross  section  of  radar  targets 

standard  deviation  of  pulse  or  video  integrated-  received  power  for 
any  given  number,  n,  independent  samples 

standard  deviation  of  pulse  or  video  integrated  received  power  for 
single  independent  samples 

"reshuffling  time",  or  time  required  for  hydrometeors  in  the  radar 
pulse  volume  to  redistribute  themselves  sufficiently  to  give  non- 
autocorrelated,  independent  data 

factor  specifying  fraction  of  radar  pulse  volume  filled  with  scatterers 

constant  of  proportionality  in  standard  deviation  versus  number  of 
independent  samples  equation 

general  symbol  for  solid  angle 


Other  S\mlx>l> 

i _ i  when  used  to  "underline"  any  given  Arabic  or  (Ireek  symbol,  represents 

the  decibel  equivalent  value  of  the  parameter.  For  example, 

JL  10  log  x. 
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